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Abstract 



A precise determination of the energy scale of jets at the Colhder Detector at 
Fermilab at the Tevatron pp collider is described. Jets are used in many analyses to 
estimate the energies of partons resulting from the underlying physics process. Sev- 
eral correction factors are developed to estimate the original parton energy from the 
observed jet energy in the calorimeter. The jet energy response is compared between 
data and Monte Carlo simulation for various physics processes, and systematic un- 
certainties on the jet energy scale are determined. For jets with transverse momenta 
above 50 GeV the jet energy scale is determined with a 3% systematic uncertainty. 
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1 Introduction 



Measurements of hard scattering processes in pp collisions often depend on the deter- 
mination of the four-momenta of quarks and gluons produced in the hard scatter. The 
measurement of these four-momenta relies on the reconstruction of hadronic jets, resulting 
from the quark or gluon fragmentation. 

At the Collider Detector at Fermilab (CDF) jets are observed as clustered energy 
depositions in the calorimeters. In this article we describe how these jets are then corrected 
to correspond to the energy of the parent parton. The precision to which this can be 
achieved determines the precision of many measurements, e.g. a 1% uncertainty on the 
energy scale of jets results in an uncertainty of 10 % on the cross section for jet production 
at transverse momenta of 500 GeV/c ^ and in a 1 GeV/c^ uncertainty on the top quark 
mass 12 • 

The original parton transverse energy can be estimated by correcting the jet for in- 
strumental effects and for radiation and fragmentation effects: 

parton _ ^^jet (j^^^^^ x C^bs — Cue + Cqoc = p^"''*"^'^ — Cue + Cooc (1) 

where is the transverse momentum of the parent parton the procedure is aimed at, 

p^^* is the transverse momentum measured in the calorimeter jet, p^f'^**'^'^ is the transverse 
momentum of the particle jet, that is, a jet corrected by all instrumental effects which 
corresponds to the sum of the momenta of the hadrons, leptons, and photons within the 
jet cone, and 

• Crj, "?7-dependent" correction, ensures homogeneous response over the entire angular 
range; 

• Cmi, "Multiple Interaction" correction, is the energy to subtract from the jet due 
to pile-up of multiple pp interactions in the same bunch crossing; 

• CAbs, "Absolute" correction, is the correction of the calorimeter response to the 
momentum of the particle jet. Particle jets can be compared directly to data from 
other experiments or theoretical predictions which include parton radiation and 
hadronization. 

• Cue and Cqoc, the "Underlying Event" and "Out-Of-Cone" corrections, correct for 
parton radiation and hadronization effects due to the finite size of the jet cone algo- 
rithm that is used. Note that these corrections are independent of the experimental 
setup, i.e. the CDF detector environment. 

These corrections and their systematic uncertainties will be described below and in 
the following sections. All the correction factors are determined as a function of the jet 
transverse momentum but they apply to all components of the four-momentum of the jet. 
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The CAbs correction is derived using a detailed Monte Carlo (MC) simulation of the 
physics processes and the detector response. The corrections Cjje and Cqoc are deter- 
mined using the PYTHIA MC generator. Thus the major task is the tuning and validation 
of the detector simulation as well as of the physics modeling used in the simulation. The 
other corrections are mostly derived directly from data but are also compared to the 
simulation. 

The following ingredients are necessary for deriving the corrections described above: 

• The energy scale for the electromagnetic calorimeter is set using electrons from the 
decay Z e~^e~ . The energy scale for the hadronic calorimeter is set to the test- 
beam scale of 50 GeV/c charged pions. Section |21 describes the CDF detector and 
the definition of the calorimeter energy scales. 

• Jets are defined using a cone algorithm either on calorimeter towers, on stable 
particles, or on partons in the MC. Different cone sizes are studied and all the 
corrections are derived for the specific cone size. The details of the algorithm used 
for defining a jet are given in Sec. El 

• Many datasets are used for either developing a correction procedure, or for validating 
or tuning the MC simulation. These are described in Sec. |3] 

• Since the simulation is used to correlate a particle jet to a calorimeter jet a detailed 
understanding of the detector simulation is needed. Therefore the simulation is 
tuned to model the response of the calorimeter to single particles by comparing 
the calorimeter energy measurement, E, to the particle momentum, p, measured in 
tracking detectors. Here, measurements based on both test beam and CDF data 
taken during Run II are used. The details of the simulation are given in Sec. El 

• The calorimeter simulation is most reliable in the central part of the calorimeters 
since the tracking coverage in the forward regions is limited. Therefore, the forward 
calorimeter jet response is calibrated with respect to the central, to flatten out 
the jet response versus the jet polar angle. This procedure also corrects for the 
lower response in poorly instrumented regions of the calorimeters. The ?7-dependent 
correction, C^, is described in Sec. (HI 

• After tuning the simulation to the individual particles response and achieving a 
jet response independent of the polar angle, calorimeter jets are corrected to a 
particle jet, i.e. they are corrected for the central calorimeter response. The absolute 
correction, CAbs, is derived from the simulation and described in Sec. [71 Since 
the correction is derived from simulation, it is also important to ensure that the 
multiplicity and momentum spectrum of particles in the data is well reproduced by 
the simulation. This is also presented in Sec. [71 
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• A further correction is made for pile-up of additional pp interactions. This pile-up 
can lead to an overestimate of the jet energy if particles produced in the additional 
interactions happen to overlap those produced in the hard scattering process. A 
correction, Cmi, for this is derived from data and described in Sec. |H1 

• The jet energy needs to be corrected for particles from the underlying event, i.e. 
interactions from spectator quarks and initial state QCD radiation, since the mea- 
surement aims at estimating the parton energy. This correction. Cue, is described 
in Sec. M 

• Since the jet cone is of finite size some particles originating from the initial parton 
may escape from the jet cone either in the fragmentation process or due to parton 
radiation. The out-of-cone energy, Cqoc, is measured in MC simulated events and 
compared to the data. This correction is described in Sec. IHl 

• Various cross-checks using different physics processes are presented to validate the 
universality of the procedure and verify the systematic uncertainties. These are 
presented in Sec. Uni 

• A summary of the systematic uncertainties is given in Sec. El They take into 
account any differences observed between the data and the simulation and possible 
systematic biases in the procedure used to determine the corrections. 

In Sec. El we present conclusions and an outlook of future possible improvements on 
this correction procedure and the systematic uncertainties. 
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2 CDF Detector 



The CDF Run II detector has been described in detail elsewhere jS]. In the following we 
use a cylindrical coordinate system with the origin at the center of the detector where the 
z axis points along the beam pipe in which Q is the polar angle, is the azimuthal angle 
and T] = — lntan(^/2) is the pseudo-rapidity. The transverse energy, Et, is defined as 
i?sin6' and the transverse momentum, pt, as psin^, where E is the energy measured by 
the calorimeter and p the momentum measured in the tracking system. The imbalance 
in transverse energy, I^t, is the magnitude of I^t with I^t = —Y^i E^^i, where Ui is a 
unit vector that points from the interaction vertex to the ith calorimeter tower in the 
transverse plane. 

Transverse momenta of charged particles {px) are measured by an eight-layer silicon 
strip detector j3] and a 96-layer drift chamber inside a 1.4 Tesla magnetic field. The 
innermost layer of the silicon detector is located on the beam pipe at a radius of 1.5 cm, 
with the outermost layer located at 28 cm. The silicon detector provides tracking in the 
pseudo-rapidity region |?7| < 2, with partial coverage up to |?7| < 2.8. Outside of the silicon 
detector, the Central Outer Tracker (COT) jS] is a 3.1 m long, open-cell drift chamber 
with an active tracking region extending radially from 41 cm to 137 cm. The COT's 
96 layers are divided into super-layers of 12 wires each that alternate between axial and 
stereo orientation. The COT provides coverage for |?7| < 1. The efficiency for finding 
charged particle tracks is close to 100% for \ri\ < 1 and falls to about 40% for |?7| ^ 2. 
The momentum resolution is (t{pt)/pt = 0.15% x pt for |?7| < 1.0 and degrades with 
increasing |?7|. 

Located outside the solenoid, a segmented sampling calorimeter is installed for the 
measurement of the electromagnetic and the hadronic energy depositions, which is de- 
scribed in detail in Sec. 12.11 The central and forward part of the calorimeter have their 
own shower profile detector positioned at the expected maximum of the lateral shower 
profile, the Central Electromagnetic Showermax (CES [B]) and the Plug Electromagnetic 
Showermax (PES [7j) detectors. Located at the inner face of the central calorimeter, the 
Central Pre- Radiator (CPR jS]) chambers use the solenoid coil as a radiator to measure 
the shower development. These three detectors are mainly used for photon and elec- 
tron identification. Drift chambers located outside the central calorimeters and detectors 
behind a 60 cm iron shield detect energy depositions from muons with \rj\ < 0.6 jH]. Addi- 
tional drift chambers and scintillation counters detect muons in the region 0.6 < {rjl < 1.0. 
Luminosity monitoring is provided by the Cherenkov Luminosity Counter (CLC) |10j . 
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2.1 Calorimeters 



The CDF calorimeter is divided into a central and a forward section. A schematic view 
is shown in Fig. [T] 




Figure 1: Elevation view of one half of the CDF detector displaying the components of 
the CDF calorimeter: CEM, CHA, WHA, PEM and PHA. 

There are a total of five calorimeter compartments: the central electromagnetic, 
CEM and central hadronic, CHA ^2], the plug electromagnetic, PEM ^3], and 
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plug hadronic, PHA jT2], and the wall hadronic WHA jT^j in the region between the 
central and the forward calorimeter a.t 6 ^ ±30° (see Fig. 

The central calorimeter covers the region |?7| < 1.1 and is divided in two halves at 
|?7| = 0. It is segmented in towers of 15° in azimuth and 0.1 in t] with lead-scintillator sam- 
pling for the electromagnetic measurements and steel-scintillator sampling for hadronic 
measurements. The material in the CEM has a depth of 18 radiation lengths. The en- 
ergy resolution for high energy electrons and photons is ^^-^ = ^^^^ © 1.5%. The CHA 
and WHA are of similar construction, with alternating layers of steel and scintillator 
and are 4.7 interaction lengths deep. Three of the WHA towers are situated behind the 
CEM / CHA and three are behind the plug calorimeter. The energy resolution of the CHA 
is = ^ ©3% and WHA is = ^ ©4% for charged pions that do not interact 
in the CEM. ^ 

The forward "plug" calorimeters cover the angular range corresponding to 1.1 < |?7| < 
3.6. They are segmented in 7.5° towers for |?7| < 2.11 and 15° for \r]\ > 2.11. The PEM 
has a depth of 23.2 radiation lengths. The energy resolution for high energy electrons and 
photons is = © 1%. The PHA has alternating layers of iron and scintillating tile, 

for a total of 6.8 interaction lengths. The energy resolution of the PHA is ^ = ^ ©5% 



for charged pions that do not interact in the PEM. 

Each calorimeter tower is read out by two photomultipliers in the CEM, CHA and 
WHA and by one photomultiplier in the PEM and PHA. The calorimeter readout elec- 
tronics was upgraded for the Run II data taking period to accommodate the 396 ns beam 
bunch spacing as well as a possible upgrade to 132 ns. The ADC integration gate for the 
charge collection is 120 ns wide and collects 94 to 98% of the signal depending on the 
calorimeter type. During test beam and the previous Run I data taking period (1992- 
1996) the time between two bunches was 3.5 /is and the integration time was 600 ns for all 
calorimeter compartments. The fractional energy loss due to the shorter ADC integration 
gate in Run II is measured in muon, electron and jet data with an uncertainty of 1.5%. 



2.2 Definition of the Energy Scale of Calorimeter Towers 

The absolute energy scale of the CEM calorimeter is set such that for fully corrected 
electrons ^3] the measured mass of the Z boson in the e~^e~ decay mode is consistent 
with the mass measured at LEP ^3] taking into account photon radiation. The ratio of 
the measured calorimeter energies and the track momenta for electron candidates, E/p, 
is used to apply additional relative calibrations for each tower to improve the resolution 
of the energy measurement. The PEM energy scale is set using Z —>■ e^e~ events with 
one electron in the CEM and one electron in the PEM. 

For the hadronic calorimeters, CHA, WHA, and PHA, the initial energy scale is defined 
by their responses to a charged pion test beam of 50 GeV/c (see Sec. Eland IT7| ITH]) 
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using pions with almost no interaction in the respective electromagnetic compartments 
CEM and PEM, respectively. The total raw energy deposited in a tower is given by the 
sum of CEM and CHA energies. For charged pions that interact in the CEM the response 
is lower (see Sees. |S]and[7j). 

2.3 Stability of the Energy Scale 

The stability of the CDF calorimeter is monitored online using various calibration meth- 
ods. The energy scale of both the electromagnetic and hadronic calorimeters in general 
decrease with time due to aging of both the scintillators and the photomultipliers. The 
online response is kept stable to better than 3% while offline a stability better than 
0.3 — 1.5% is achieved. The following methods are used to obtain this stability: 

• In the CEM, E/p oi electrons with transverse energies, Et > 8 GeV are used to 
monitor the time dependence. The energy scale decreases by 3% every 6 months, 
and is corrected accordingly. 

• In the CHA and WHA the energy scale is monitored using three independent meth- 
ods: a laser system ^2]; muons from J/ip ^ fi'^fi~ decays and minimum bias data 
(these datasets are defined in Sec. |^. The test beam energy scale has been main- 
tained since 1987 using Cs^^^ source calibration runs. The CHA response decreases 
by about 1% and the WHA response by about 3% per year which is corrected by 
adjusting the calibration. 

• The PEM and PHA calorimeters are monitored using a laser system jSl and a 
radioactive source calibration using Co^°. The laser is only sensitive to aging of 
the photomultipliers while the source is sensitive to both the photomultiplier and 
scintillator aging. The plug calorimeter scale decreases by up to 2-10% per year, for 
|?7| = 1.2 — 3.6. The largest decrease is observed in the region closest to the beam 
pipe. This decrease is calibrated accordingly. 

The calibration stability of the electromagnetic scales is verified using the time de- 
pendence of the reconstructed invariant Z boson mass. Figure |21 shows the mean of the 
Z e~^e~ mass distribution between 86 and 98 GeV/c^ as a function of the run num- 
ber. The range of run numbers corresponds to the data taking period from April 2002 
until September 2004. The mean energy deposited by muons from W — » /iz/^ candidate 
events is shown in Fig. 01 for muons with |?7| < 1 to verify the stability of the CHA and 
WHA energy scales. For the PHA both muons and the jet response are used to verify the 
stability. 
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Figure 2: Mean invariant mass of Z ^ e^e candidates, (M(ee)), versus run number 
for events with 86 < M(ee) < 98 GeV/c^. Shown are the values for events with both 
electrons in the central calorimeter (full circles) and for events with one electron in the 
central and one in the plug calorimeter (open circles). The dashed lines indicate a ±0.3% 
uncertainty around 91.1 GeV/c^. 



From Figs. |21and|nian uncertainty on the stability of the CEM of 0.3% and for the 
CHA of 1.5% is assigned. For a jet, which deposits typically 70% of the energy in the 
CEM and 30% in the CHA, it is thus 0.5%. 
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Figure 3: Mean energy observed in the CHA/WHA for CMUP and CMX nmons with 
Pt >20 GeV/c from W — > candidate events versus run number. The CMUP muons 
are confined to |r;| < 0.6 and thus only sensitive to the central part of the CHA. The 
CMX muons cover the region 0.6 < |?7| < 1.0 and probe the outer part of the CHA plus 
the innermost part of the WHA. The dashed lines indicate a 1.5% uncertainty. 
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3 Jet Clustering Algorithm 



The energy of a jet is calculated from the energy deposited in the calorimeter towers using 
different types of clustering algorithms. For this study, jets are clustered using a cone 
algorithm with a fixed cone size in which the center of the jet is defined as (77-^^*, 0-^'^*) and 
the size of the jet cone as i? = ^(77*"""='' - r/J<=*)2 + [(ptower _ ^j^ety ^ 9.4, 0.7, or 1.0. The 
jet corrections and uncertainties are estimated for these three cone sizes. 

3.1 Calorimeter Jets 

The jet clustering algorithm groups calorimeter towers with E^i > 1 GeV into jets. 
Eti = Ei sin^j is the transverse energy of a tower with respect to the 2;-position of the pp 
interaction, and the energy Ei is the sum of the energies measured in the electromagnetic 
and hadronic compartments of that tower. 

Firstly "seed towers" are defined in order of decreasing Exi- For each seed tower the 
towers within a radius of size R with respect to its position are used to build "clusters" . 
Once we have an initial list of clusters, the cluster transverse energy and the location of 
the cluster is calculated using the definitions: 

Ntow 

E^* = ETi (2) 

Ntow TP J. 

= E 9# (3) 

1=0 

V" = E ^ (4) 

where Nt^w is the number of towers inside the radius R with E^ > 1 GeV. 

This procedure is repeated iteratively, a new list of towers around the new center is 
determined. The jet Ej- and direction are recalculated until the list of towers assigned to 
the clusters is stable, that is, when the geometrical center of the tower correspond to the 
cluster centroid. Overlapping jets are merged if they overlap by more than 50%. If the 
overlap is smaller than 50%, each tower in the overlap region is assigned to the nearest 
jet. 

The final jet energy and momentum coordinates are computed from the final list of 
towers: 

Ntow 

Ej,t =Y.Ei (5) 

i=0 
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Ntow 

p>:'= J2 ^.sin(ft,)cos(</.,) (6) 

4 = 

Ntow 

P^y'= E ^.sin(^,)sin(0,) (7) 

j=0 

Ntow 

p^;'= J2 ^.cos(^,) (8) 

■t=0 



= \/{p'^y + ipiy (9) 

0jet = tan^ (10) 

Px 

jet 

sme,et = ^ (11) 

Erjet = E jet sin 6 jet (12) 
In general, jets with Et < 3 GeV are not used in physics analyses at CDF. 

3.2 Particle Jets in Monte Carlo 

In Monte Carlo (MC) simulation, particle jets are obtained using the same jet clustering 
algorithm on stable final state particles^, i.e. the sums in Eqns. l21IHlgo over the stable 
particles instead of the towers. A particle jet includes any particles produced in the inter- 
action, i.e. particles from the hard scattering process and also those from the underlying 
event. 



^In PYTHIA and HERWIG, these are particles which have a Ufetimc long enough to not decay within the 
CDF detector volume. 
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4 Data and Monte Carlo Samples 



In this section we describe the data and Monte Carlo (MC) simulation samples used 
throughout this paper for studying and determining the jet response. 
The following samples are used: 

• Test beam data: The response of all calorimeters was first measured in a test 
beam. Test beam for the CEM, CHA, and WHA were taken in 1985 and 1988 
[THl [T7j with a momentum range 5 — 180 GeV/ c for electrons and 7 — 220 GeV/c for 
charged pions. The test beam for the plug calorimeters was taken in 1996 ^H] with 
a momentum range 5.3 — 181 GeV/c for electrons and 8.6 — 231 GeV/c for charged 
pions. 

• Minimum bias: This sample is collected requiring at least one pp interaction. It is 
triggered by activity in the luminosity counters, i.e. the trigger requires coincident 
hits in both the east and west CLC detectors. It is used for studying multiple 
interactions and the underlying event (see Sees. IHlandlH)). Tracks from this sample 
are also used for tuning the simulation (see Sec. Ej). 

• Single track: A special trigger was designed to take data with a high momentum 
track at pr thresholds of 3, 7 and 10 GeV/c. These data are used to tune the 
central calorimeter response at larger p^ (see Sec. . The tracks in this dataset are 
confined to the region in the COT where particles can traverse all available layers, 
that is \ri\ < 1. 

• Jet: There are four samples triggered on at least one jet with transverse energy 
E^T^^ > 20, 50, 70 and 100 GeV referred to as jet-20, jet-50, jet-70 and jet-100, 
respectively. They are primarily used to calibrate the response as function of pseudo- 
rapidity (see Sec. EJ. Apart from the jet-100 sample the data are taken with a fixed 
"prescale" depending on the pt threshold. A prescale of e.g. 10 means that only 1 
in 10 events passing the trigger requirements is recorded. There is another sample 
triggered on at least one trigger tower with Et >5 GeV referred to as single-tower-5. 

• 7-jet: This sample is triggered on an isolated electromagnetic cluster with Ej, > 25 
GeV. It is used to cross-check many aspects of the calibration and to determine 
the systematic uncertainty on the out-of-cone correction (see Sec. The photon 
selection requirements are described in detail elsewhere |29j . 

• W — > eue and W — > fiu^ : These data are taken by an inclusive electron (e) 
or muon (/i) trigger with E^ > 18 GeV and pj. > 18 GeV/c, respectively. The 
electron and muon identification is described in detail elsewhere jSHl- The W — > efg 
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( W — > fiu^) sample is selected by requiring one electron (muon) with > 25 GeV 
(pj. > 25 GeV/c) and missing transverse energy I^x > 25 GeV. 

• Z — > e~^e~ and Z — >■ These data are taken with the same trigger to the 
previous W samples. They are selected by requiring two electrons (muons) with 
E!f. > 18 GeV {pj. > 20 GeV/c) and requiring the invariant mass of the electrons 
and muons to be between 76 and 106 GeV/c^. 

• J I'll) — > e+e" and J/i/^ — > The J/ip ^ e^e~ data are taken with a 
dielectron trigger with > 2 GeV for each electron. The J/ip fi~^fi~ data are 
taken with a dimuon trigger with pj, > 1.5 GeV/c for each muon. The dilepton 
invariant mass is required to be between 2.5 and 3.5 GeV/c^ for the e~^e~ decay and 
between 3.0 and 3.2 GeV/c^ for the /i"*"//" decay. 

Corresponding MC samples are used for all processes. In all cases, Monte Carlo 
samples are generated using both PYTHIA 6.216 [20] and HERWIG 6.505 [H] with CTEQ5L 
[22] parton distribution functions. 

The following MC samples are generated: 

• Minimum bias and jet: The inclusive 2^2 parton processes {pp qq + X, 
pp gq + X, pp ^ gg + X) are generated at different thresholds for the transverse 
momentum of the outgoing partons, pr > 0, 10, 18, 40, 60, 90, 150, 200, 300, 400, 
500, 600 GeV/c. The sample with > is referred to as minimum bias MC. 

• Single-particle: Single charged particles are generated with a mixture of 60% tt^, 
30% and 10% p and p. 

• 7-jet: The processes qg ^ q'-f + X and qq ^ ^g + X are generated at different 
thresholds of the transverse momentum of the outgoing partons, px > 12, 20 and 
40 GeV/c. 

• W and Z: W and Drell-Yan production is generated in both electron and muon 
decay channels as described in pUj . 

• J /ijj — > e~^e~ and J/i/^ — > This sample is generated using BGEN 

• ti — > WbWb: Top pair production is generated using HERWIG. More details can 
be found in 

• W^+jets: The ALPGEN 1 . 3 generator interfaced with HERWIG is used to simulate 
these type of events [5U] . 
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The matrix element of the hard scattering process, the underlying event, the higher 
order QCD correction (i.e. the parton showering process), and the fragmentation are all 
handled by PYTHIA and HERWIG. The matrix element is well known at leading order in 
QCD. However, the other three components are not fully calculable and thus modeled 
using an empirical approach. Note that PYTHIA and HERWIG use different models, which 
is helpful to study and understand the uncertainties associated with these effects. The 
modeling of the fragmentation and parton showering is mostly tuned to e+e~ data. For 
modeling the underlying event we use a tuning optimized to describe CDF data from 
Run-I, which we refer to as "PYTHIA Tune A" j21j. For HERWIG the default parameters 
for the underlying event are used. The recent model for the underlying event, JIMMY |22], 
has not been used in this analysis. 

After generation, the samples are processed through the CDF Run II detector si- 
mulation, which is described in detail in Sec. El They are then processed through the 
standard reconstruction program as the data. The details of particle and calorimeter jet 
reconstruction are in Sec. El 

Throughout this article the data are primarily compared to PYTHIA MC samples but 
any differences between PYTHIA and HERWIG are discussed. Detailed comparisons between 
HERWIG and PYTHIA are presented in Sec. [TUl 
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5 Calorimeter Simulation 



In this section we describe the simulation of the CDF calorimeter. A parameterized 
shower simulation is used with the simulation parameters tuned to describe the observed 
calorimeter response of single particles. The response of single isolated charged particles is 
measured using test beam data as well as minimum bias and single track trigger data from 
Run II. Similarly, the tuning of electromagnetic showers is based on electrons observed 
in J /if) — i> e"'"e~ and Z — > e+e~ decays. In the following, we first describe the shower 
parameterization, then the measurement of the calorimeter response is discussed, and 
finally the systematic uncertainty is assigned based on how well the MC simulation models 
the data. 

5.1 Simulation of Electromagnetic and Hadronic Showers 

The CDF detector simulation uses GEANT ^2Z| to track generated particles through the 
CDF detector and to simulate secondary physical processes such as energy loss, multiple 
scattering, and inelastic interactions. After the first inelastic interaction in the calorime- 
ter, the particles are passed to the program GFLASH j2H|; a fast simulation of electromag- 
netic and hadronic particle shower. GFLASH generates particle showers shapes within the 
calorimeter and computes the energy deposited in the calorimeter sensitive volumes, using 
parameterizations of the longitudinal and lateral shower profiles. The parameterizations 
of electromagnetic and hadronic showers are described in detail in Ref. [2H] and are briefly 
outlined in this section. 

5.1.1 Procedure 

The simulation of electromagnetic and hadronic showers involves two steps. First, GFLASH 
calculates the spatial distribution of energy, Edp-, deposited by a shower within the calorime- 
ter volume: 



where the longitudinal energy profile, L{z), depends on the shower depth z, and the lateral 
shower profile, T(r), depends on the radial distance, r, from the trajectory of the particle. 
This parameterization takes dependencies on the incident particle energy and shower 
fluctuations into account, and further considers the repetitive sampling structure of the 
detector volume. Second, the fraction of the deposited energy which is visible to the active 
medium, Emsif), is determined. Emsif) is computed depending on the relative sampling 
fractions of minimum ionizing particles, electromagnetic and hadronic particles. The 
relative sampling fractions for electrons, Se/ Smip, and for hadrons, Shad/Smip, compared 




(13) 
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to the sampling fraction of a particle that does not interact inelastically in the calorimeter, 
Smip, are two tunable parameters. 

5.1.2 Longitudinal shower profile 

Electromagnetic Showers: For the simulation of the longitudinal profile of electro- 
magnetic showers GFLASH assumes that they follow the F-distribution: 



Lem{z) = —— — X = PemZ{Xo) (14) 

1 (aem) 

where z is the shower depth, measured in units of radiation lengths, Xq. The parameters 
aem and /3em are Gaussian-distributed, and parameterized as a function of incident particle 
energy. The longitudinal profile of each shower is simulated by choosing a correlated 
{aem, Pern) pair at random. 

Hadronic Showers: GFLASH distinguishes three classes of hadronic showers: 

• Purely hadronic showers (h) whose propagation scales with the absorption length 
Ao. 

• Showers where a 7r° is produced in the first inelastic interaction (/). Their propa- 
gation scales with radiation length Xq. 

• Shower profiles from tt'^'s produced at a later stage of the hadronic showering process 

(0- 

The hadronic showers in GFLASH are calculated as a superposition of those three shower 
classes: 



with 



dEdp = fdp Einc [chLh{xh)dxh + CfLf{xf)dxf + ci Li{xi) dxi] (15) 



Ci^{xh) = '^h^l^ ^^^p^^^Xq) (16) 

'^f~^ -Xf 

= Xf = Pfz{Xo) (17) 

= ^Vrr- ^' = A^(Ao) (18) 



(19) 
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The coefficients Ch, Cf and q are the relative fractions of the three contributions nor- 
mahzed such that Ch + cj + q = 1. The factor f^p is the fraction of deposited energy 
with respect to the energy of the incident particle, and takes the intrinsic losses during 
the hadronic shower development into account. 

The relative probabilities of the three classes Ch, Cf and q depend on the incident 
energy and are correlated through: 

Ch{E) = 1 - fAE), Cf{E) = fAE)a - fAE)), ci{E) = fAE)fAE) (20) 

where /ttO {E) is the probability that a hadronic shower contains any and {E) is the 
probability that a 7r° is produced "late", i.e. not in the ffist interaction. 

In total, the longitudinal hadronic shower development depends on 18 partially corre- 
lated parameters: the mean and a values of the three different a- and /3- values for each 
of the three shower types, and the fractions fdp, fnO and flo- For the electromagnetic 
shower development there are only 4 parameters: the mean and a values of aem and Pem- 



5.1.3 Lateral shower profile 

The parameterization for the lateral energy proffie of electromagnetic and hadronic show- 
ers is taken to be 

(r^ + RioV 

-R50 is given in the units of Moliere radius Rm for electromagnetic showers, and absorption 
length, Aq, for hadronic showers. -R50 is an approximate log- normal distribution, with a 
mean value {R^q{E,z)) and variance a^^g parameterized as a function of the incident 
particle energy, E, and the shower depth, z, 

{RAE^z)) = [Ri + iR2-R3\ogE)zr n=l,2 (22) 
anAE,z) = [(^1 + (82 - SslogE)z{RoiE,z))f (23) 

The ^-evolution of the lateral spreading is linear for hadronic showers (n = 1) and 
quadratic for electromagnetic showers (n = 2). The electromagnetic and hadronic proffies 
are determined by their own set of adjustable Ri and Si values, thus giving a total of 14 
parameters. 



5.1.4 Tuning to CDF Data 

The GFLASH longitudinal and lateral hadronic shower parameters were tuned using single, 
isolated tracks selected from minimum bias data (0.5 < p < 2.5 GeV/c in the central 
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and 0.5 < p < 5 GeV/c for the plug calorimeter) and test beam data {7 < p < 220 
GeV/c). The electromagnetic showering parameters were tuned using test beam data 
and compared in situ using Z — > e^e~ events data and simulation. Further adjustments 
in addition to the above parameters are made to the energy deposited at the boundary 
between central and plug calorimeter. 

Given the limited availability of CDF isolated single track and test beam data with 
current statistics, not all of the 34 parameters described above were tuned to CDF data. 
Furthermore, the current tuning of various parameters based on CDF data is restricted 
to relatively low particle momenta, below 2.5 GeV/c. For the remaining parameters, we 
use the default setting from the HI collaboration jSHj- At this stage the parameters that 
are tuned to CDF data are: 

• Oh and ai for both the central and plug calorimeters. 

• /?/ for both the central and plug calorimeters. 

• fdp and J'tjO for both the central and plug calorimeters. 

• four of the parameters that characterize the lateral profile for the central calorimeter: 
Ri, i?2, Si and 5*2. 

• the relative sampling fractions Se/Smip and Shad/Smip for the plug calorimeters. 

Tuning these parameters gives a good description of the single track data as will be 
seen in the following sections. However, since the tuning was done on a limited dataset 
there are some discrepancies with the latest single track data samples that extend up 
to momenta of 20 GeV/c. The newly collected data that are presented in the following 
sections will be used in the near future to further tune the simulation parameters. This 
is particularly important for the plug calorimeters. 

5.2 Calorimeter Response to Hadronic Particles 

Charged and stable neutral hadrons carry approximately 70% of the jet energy and there- 
fore a good description of their response is pivotal for a good simulation of the jet en- 
ergy measurement. The hadronic shower development can best be studied using isolated 
charged particles. 

In the rapidity range covered by the central calorimeters, CDF has an excellent track- 
ing coverage and efficiency. The measurement of the response in the plug calorimeters is 
more difficult since there is no track trigger available and the tracking efficiency is limited. 
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5.2.1 Procedure 



Tracks are selected and extrapolated to the impact point at the position of the CES or 
PES detectors, taking the magnetic field into account. If a track is extrapolated to a 
calorimeter tower, this tower is considered to be the target tower. Tracks are required to 
point to the inner 0.9 x 0.9 contour of the target tower to stay away from the tower edges, 
otherwise they are not considered. 

From the eight towers surrounding the target tower, the measurement of the CEM 
energy uses a 2 x 2 subset containing the target tower and those three adjacent towers 
closest to the track impact point. For the CHA the shower spread is typically larger and 
all eight adjacent towers are considered in the signal energy definition. A sketch of the 
target tower and signal region definitions is shown in Fig. |3]for the CEM and the CHA. 



EM: 2x2 HAD: 3x3 




Figure 4: Illustration of the target tower for the electromagnetic (EM) and the hadronic 
(HAD) sections used in the tuning of charged hadrons. The yellow (lightest shading) 
region is the target tower. The "x" marks the impact point of the track. The yellow and 
cyan (light and medium shading) region is the signal region and the dark blue (darkest 
shading) illustrates the background region. The horizontal axis represents the t] direction 
and the vertical represents the direction. 

The signal tower region may contain additional particles overlapping with the direction 
of the primary track that is being analyzed. This background is reduced by requiring no 
additional tracks within a 7 x 7 block of towers around the target tower. In addition, no 
energy depositions above 1 GeV in the CES detector are allowed within the 7x7 blocks 
of towers except within AR = ^ {r]cES - VtrackY + [(f^CES - (ptrackf < 0.03 around the 
extrapolated track position, where (pcES {4>track) and ijces {vtrack) denote the azimuthal 
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angle and the pseudo-rapidity of the CES cluster (track). The first requirement removes 
charged particles and the second requirement mostly removes photons. 

Even after these isolation requirements, some overlap background remains. The back- 
ground is estimated using the energy deposited in the same r] range but within the towers 
along the edge in of a 5 x 5 tower group around the target tower, using those edge 
which is more distant from the track impact point, see Fig. |3] The energy measured in 
those background regions is then scaled to the area covered by the signal region and is 
subtracted from the energy measured in the signal regions. For the single track MC the 
effect of the background subtraction is negligible. 

5.2.2 Central Calorimeter 

The comparison of the mean value of the E/p distribution, (E/p), between single track 
data and MC is shown in Fig. El The data are shown before and after background 
subtraction. The background contributes significantly for p < 3 GeV/c but is negligible 
at higher p. For particle momenta increasing from 0.5 to 20 GeV/c, the fraction of energy 
deposited in the CEM drops from 40% to 25%, whereas the fraction in the CHA rises 
from 20% to 55%. In total {E/p) rises from about 0.5 ai p = 0.5 GeV/c to about 0.8 
at p > 5 GeV/c. For p > 5 GeV/c the response is almost independent of p. In general, 
the mean {E/p) agrees well between the data and the simulation separately for the CEM 
and the CHA energies and for their sum. The deviations observed at very low momenta 
and around 4 GeV/c are probably related to differences in the particle spectrum between 
data and MC caused by momentum cutoffs in the MC and trigger thresholds in the data. 

Beyond 20 GeV/c there are not enough isolated tracks to verify the simulation and 
test beam data are used for studies of higher momenta. The simulated CEM, CHA, and 
CEM-I-CHA energy is compared with the test beam response for 57 GeV/c pions in Fig. 
|B] For this study the charged pion sample is divided into two categories: 

• pions which do not interact in the CEM are selected by requiring the CEM energy 
to be less than 500 MeV ^. The CHA energy of these pions is shown in the upper 
left display in Fig. 13 This category is also used to set the absolute energy scale of 
the hadronic calorimeters (see Sec. [7j). 

• pions which interact inelastically in the CEM: the CEM, CHA and total energy of 
this category is displayed in the other three plots in Fig. [HI 

The data are reasonably well described by the simulation. 

Fig. [7| shows a summary plot of the absolute total response {E/p) for minimum bias, 
single track trigger and test beam data together with the Monte Carlo expectation. 

^The mean CEM response for non-interacting particles is about 300 MeV. 
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Figure 5: Fractional energy observed in the central calorimeter as a function of incident 
particle momenta. The top row shows {Ecem/p), {Echa/p) and {{Eqem + Echa) Ip) for 
data signal (triangles) and background (histogram) and for single track MC simulation 
(open circles). The bottom row shows the same distributions for data after background 
subtraction (full circles) and MC simulation (open circles). 



Note, that in test beam the particles were shot at the center of the tower while for 
the single track data the inner 81% of the towers are used. This difference in selection 
is estimated to cause about a 5% higher response for the test beam data. In addition, 
for the test beam only one central calorimeter tower (covering the region 0.3 < |?7| < 0.4) 
was used. Therefore, the data cannot be compared directly between test beam and in 
situ measurements. However, in each case the results are consistent between data and 
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Figure 6: Top left: Energy observed in CHA for particles that do not interact in the CEM. 
Total (top right), CEM (bottom right) and CHA (bottom left) energy for charged pions 
with p = 57 GeV/c. The test beam data (points) are compared to the CDF simulation 
(solid line). 



simulation. 
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Figure 7: (E/p) for the central calorimeter versus particle momentum p. The top plot 
shows the in situ minimum bias (circles) and single track data (squares) compared to 
the simulation (open symbols). The bottom plot compares the test beam data (closed 
triangles) to the simulation (open triangles). The error bars on the minimum bias and 
single track data are statistical. For the test beam data we display an uncertainty of 3.5% 
due to the systematic uncertainty on the test beam momentum scale which is correlated 
between all the data points. For the test beam comparison the MC points have been 
shifted slightly in p-value to allow a better visual comparison. 
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5.2.3 Plug Calorimeter 

In the plug calorimeter the same procedure is used as in the central calorimeters. However, 
due to the smaller tower size a larger fraction of the incident particle momentum is 
measured by neighboring towers. Thus, for the (E/p) analysis, two adjacent towers in 
azimuth are treated as one target tower. 

Some difficulties arise in the plug calorimeter which are not present in the central 
calorimeter: 

• The tracking efficiency is rather low in the forward region of the detector, decreasing 
from 70% at |?7| ~ 1.2 to 30% at |?7| = 2.2. Thus the background rejection is less 
efficient in the plug than in the central region where the efficiency is nearly 100%. 
Since the tracking efficiency in the very forward region |?7| > 2.2 is even lower this 
method cannot be used to evaluate the response. 

• The triggers dedicated to collect high momentum tracks are limited to the COT 
coverage |?7| < 1.0. 

• The momentum resolution for tracks is poorer than in the central region because 
the track reconstruction is mostly based on the Silicon Vertex Detector. This re- 
sults in a less accurate determination of the reference scale, p, particularly in the 
high momentum region. However, this problem can be minimized by using tracks 
with partial coverage by the COT. Tracks with combined COT and silicon hits are 
available up to |?7| < 1.8. 

A measurement of the single particle response in the plug using combined COT and sil- 
icon tracks is shown in Fig.|Hlfor single track trigger data and simulation. The background 
is larger than in the central calorimeter as expected. After background subtraction the 
data and simulation agree well at low momenta but deviations of up to 13% are observed 
between 5 and 10 GeV/c. Note that for the tuning relevant for this paper the data at 
the medium momentum range were not available. Using the data shown here the plug 
simulation will be improved in the near future. 

Figure El shows {E/p) versus p for the plug calorimeter for minimum bias data and test 
beam data compared with the corresponding simulations. At low momenta the response 
is about 60%, increasing to nearly 100% at high momenta. 

Due to the discrepancies described above the plug calorimeter simulation is not used 
for determining the absolute jet energy scale. An additional calibration is made where 
the plug calorimeter response is calibrated with respect to the response of the central 
calorimeter using dijet events as will be described in Sec. IHl Therefore, no systematic 
uncertainty is associated directly with the {E/p) measurements for the plug calorimeter. 
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Figure 8: (E/p) observed in the plug calorimeter as a function of incident particle 
momenta. The top rows shows Epem/p, Epha/p and {Epem + Epha)/p for data signal 
(closed squares) and background (solid line) and for MC simulation signal (open squares) 
and background (dashed line). The bottom row shows the same distributions for data 
after background subtraction (full circles) and MC simulation (open circles). The tuning 
of the calorimeter was done using tracks up to 5 GeV/c. 



5.3 Calorimeter Response to Electromagnetic Particles 

On average about 30% of the particles in a jet are neutral pions which mostly decay into 
two photons: vr*^ — > 77. Therefore a good understanding of the calorimeter response to 
electromagnetic particles is important. The electromagnetic single particle response is 
studied similarly to the hadronic response, using the track momentum as reference and 
EM energies measured in the target tower and the two towers adjacent in 77 around the 
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Figure 9: (E/p) versus p for the plug calorimeters. The upper plot shows a comparison 
of single track and minimum bias data (closed squares) and MC (open squares). The 
bottom plot shows the comparison of test beam data (closed triangles) and single track 
MC (open triangles). The tuning of the calorimeter was done using tracks up to 5 GeV/c. 
For the test beam comparison the MC points have been shifted slightly in p-value to allow 
a better visual comparison. 
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track impact point. 

A comparison of the (E/p) response for electrons and positrons is shown in Fig. ITUlin 
W e^z/g and J/ip e^e~ events. Overlaid is the simulated response from correspond- 
ing MC samples. A certain momentum dependence of {E/p) arises from the W and J/il) 
selection cuts and due to final state radiation of photons. The simulation reproduces the 
data to better than 1% accuracy. 
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Figure 10: {E/p) versus p for electrons and positrons from J/ip —>■ e~^e~ and W e^Ue 
data (closed triangles and circles) and MC samples (open triangles and circles) samples. 
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5.4 Uncertainties 

Figure ^2 is used to estimate the systematic uncertainty on the central calorimeter energy 
determined from the difference between data and simulation for charged hadrons: 

We have derived the following momentum dependent estimates for the average differ- 
ences: 

• 1.5% for p<12 GeV/c 

• 2.5% for 12 < p < 20 GeV/c 

• 3.5% for p > 20 GeV/c 
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Figure 11: Difference between (E/p) of data and MC versus particle momentum p. 
Sfiown is tlie difference for minimum bias Data and MC (squares), single track data and 
MC (circles) and test beam data and MC (triangles). The dashed lines indicate the quoted 
systematic uncertainties. The error bars on the minimum bias and single track data are 
statistical. For the test beam data we do not display any uncertainty. 



The uncertainties at p < 20 GeV/ c directly reflect the limited performance of the calorime- 
ter simulation as well as limited available single isolated track statistics at medium mo- 
menta. The uncertainties at p > 20 GeV/c are due to uncertainties in the test beam 
momentum scale (2%), the shorter integration time in the CDF detector readout com- 
pared to the test beam measurement (1.5%) (see Ref. [121 CZj and Sec. Adding those 
uncertainties linearly we obtain the systematic uncertainty for p > 20 GeV/c of 3.5%. 

The {E/p) measurements presented so far are only sensitive to the inner 81% of the 
tower. In particular the instrumentation between the tower 0-boundaries is limited, and 
the exact modeling of this region in the simulation is difficult. Figure IT^ shows {E/p) 
versus relative {<prei)i that is the azimuthal angle of the track impact point with respect 
to the target tower center, normalized to the of the tower edges such that (prei = 
represents the tower center and (prei = ±1 the tower boundaries. Shown are data and 
simulation for p = 3 — 5 GeV/c and p = 12 — 16 GeV/c. 

Data and MC differ by up to 10% near the 0-boundaries (If/'r-eil > 0.9). Similar 
discrepancies are seen at the rj boundaries of the towers. These 10% differences are taken 
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as systematic uncertainty. Since these boundaries in rj and correspond to 19% of the 
total tower area, the 10% uncertainty at the boundaries translates into a 1.9% systematic 
uncertainty on the overall particle response to charged hadrons. 




Single track data o Single track MC 



▼ MC - Data 



Figure 12: (E/p) vs. (prei for particles with momenta between 3 and 5 GeV/c (left) and 
for particles with momenta between 12 and 16 GeV/c (right). Top: data (full points) and 
MC (open squares). Bottom: difference between data and MC. 
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Figure ^1 shows the difference in (E/p) between data and simulation for electrons. 
The data and the simulation agree to within 1% which is taken as systematic uncertainty. 
This measurement only uses tracks pointing to the inner 84% of a tower in azimuth, (prei, 
as a consequence of the electron selection that involves a CES energy cluster fiducial cut. 
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Figure 13: Difference between {E/p) in data and simulation versus p for electrons from 
J/ip —>■ e^e~ and W — e^u^ (closed squares). The individual distributions for data and 
simulation are shown in Fig. The dashed line indicates the systematic uncertainty on 
the electromagnetic energy scale. 



The remaining 16% at the boundaries in are studied using Z — s> e^e~ events where 
one electron is identified using the standard criteria jHO] and the second electron is just 
taken to be the highest momentum track with 35 < p < 55 GeV/c in the angular range 
covered by the central calorimeter. There is no requirement on the energy measured in 
the calorimeter and thus the measurement is unbiased toward the calorimeter response. 
The invariant mass between the electron and the candidate track is required to be within 
10 GeV/c^ of the Z boson mass. For these track based electrons candidates, {E/p) is 
shown in Fig. El as function of (prei- 

Figure im also shows the difference between data and simulation in {E/p) in Z ^ e+e~ 
events as a function of relative </>. A discrepancy of about 10% is observed between 
data and simulation at \(f)rei\ > 0.84. This region corresponds to 16% of the tower area, 
thus causing an overall 1.6% difference between data and simulation for the response of 
electromagnetic particles. 

A summary of the systematic uncertainties is given in Tabled These uncertainties are 
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Figure 14: Top: (E/p) versus (prei for electrons in Z — > e~^e~ events. (E/p) is shown as 
a function of relative 0. Shown are data (closed points) and simulation (open circles). 
Bottom: The difference in {E/p) between data and simulation for electrons. 



uncorrelated and added in quadrature. The overall systematic uncertainty on the energy 
of charged hadrons is thus 2% for p < 12 GeV/c, 3% for 12 < p < 20 GeV/c and 4% for 
p > 20 GeV/c^. The total systematic uncertainty on the energy of electromagnetically 
showering particles is 1.7% for the entire momentum range studied. 
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p (GeV/c) 


0-12 12-20 >20 


(E/p) response to hadrons 




Total tower (%) 

Near tower and ?7-boundaries (%) 


1.5 2.5 3.5 
1.9 1.9 1.9 


Total for hadrons(%) 


2.5 3.0 4.0 


(E/p) response to EM particles 




Total tower (%) 

Near tower 0-boundary (%) 


1.0 1.0 1.0 
1.6 1.6 1.6 


Total for EM particles(%) 


1.7 1.7 1.7 



Table 1: Summary of the relative uncertainties due to the modeling of the simulation of 
charged hadron showers and electromagnetic (EM) showers. All the uncertainties refer 
to the effect on the entire tower energy (see text). The total uncertainty is the quadratic 
sum of the individual uncertainties weighted with the area covered by the individual 
measurements. 
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6 Ty-dependent Corrections 



Even after the calorimeter energy calibrations described in Sec. |21 the response of the 
CDF calorimeter is not uniform in pseudo-rapidity. The dependencies on t] arise from 
the separation of calorimeter components at 77 = where the two halves of the central 
calorimeter join and at |?7| ~ 1.1 where the plug and central calorimeters join. The dif- 
ferent responses of the plug and central calorimeters also cause a dependence on 77. The 
?7-dependent corrections are introduced to flatten the 77 dependence of the calorimeter re- 
sponse. The method implicitly also includes corrections for both the transverse spreading 
of calorimeter showers outside the jet cone and any rj dependence of gluon radiation and 
multiple parton interactions. 



6.1 Correction Procedure 

The ?7-dependent corrections are obtained using the "dijet balancing method". They are 
determined based on the assumption that the two jets in dijet events should be balanced 
in pt in absence of hard QCD radiation. To determine the corrections, we define a jet 
with 0.2 < |?7| < 0.6 as a "trigger jet" and define the other jet as a "probe jet". When 
both jets are in the region of 0.2 < \rijet\ < 0.6, the trigger and probe jets are assigned 
randomly. 

The Pt balancing fraction is then formed: 



A probe triqqer 

Apr ^ Pt - Vt 



where p^*^^*^^ and pj^°^^ are the transverse momenta of the trigger and probe jet, re- 
spectively. The correction factor required to correct the probe jet can then be inferred 
through 

Pdi'jet = 2 _ l^f\^ ^'^^^ 

Note, that (3dijet is mathematically equal to pj^°^'^/p^^^^^^^. However, by inferring Pdijet 
from ff, we reduce the sensitivity of our measurement to non-Gaussian tails since the 
distribution is in good approximation a Gaussian distribution unlike the distribution of 

probe I triqqer 

Pt IPt ■ 

The ?7-dependent corrections are defined as 1/ Pdijet and they are determined separately 
for data and MC and for different p^^ bins. For corrections obtained from data the samples 
single-tower-5, jet-20, jet-50, jet-70, and jet-100 are used. The corrections for the MC are 
obtained from the jet samples generated with the PYTHIA MC. 

Several cuts are placed to reduce the effects of QCD radiation: 



33 



• The jets are required to be back to back in the r — cj) plane, i.e. the difference of 
their azimuthal angle, A^Qetp^^^g, jet^^j^^g^), should be larger than 2.7 radians. 

• If a 3rd jet is present in the event, the pt of this 3rd jet should be less than 7 GeV/c 
for minimum bias data, less than 8 GeV/c for the samples triggered on 20 GeV 
or 50 GeV jets, and less than 10 GeV/c for the samples triggered on 70 GeV or 
100 GeV jets. 

• The average p^j^^ = {p^^^ +p^^*^)/2 of the two jets is required to be 5 GeV/c higher 
than the trigger threshold of the respective sample. 



• The significance of missing is defined as I^t/ \/Y,iET,i, where the sum extends 



to all the calorimeter towers. It is required to be less than 2 + 0.018 x p^'^ ■'^ 
for > 55 GeV/c and less than 3 V GeV otherwise. 

In Fig. ^1 the dijet balance (3dijet is shown for data, HERWIG and PYTHIA MC samples 
for a jet cone size of Rjet = 0.4 and four transverse momentum regions 25 < p^"^ < 55 
GeV/c, 55 < p^"^ < 75 GeV/c, 75 < p^"^ < 105 GeV/c and p^"^ > 105 GeV/c. Figures 
ITBl and IT7I show the equivalent plots for Rj^t = 0.7 and Rjet = 1-0, respectively. The lines 
in Fig. IT51IT7I show the interpolation between the individual measurements, and inverse 
of these functions is taken to be the ?7-dependent correction factor. 

It is seen that Pdijet ~ 1 in the region where the trigger jet is selected, 0.2 < |?7| < 0.6, 
for both the data and the simulation. The dips at |?7| ~ and ±1.1 are due to the gap 
between calorimeters in these regions, resulting in a lower average response. In the plug 
region, |?7| > 1.2, the calorimeter response is higher than in the central region by about 
10% at low p^"^ and 5% at high 

Both PYTHIA and HERWIG reproduce the data well up to \ri\ = 1.4 at all At larger 
|?7| and for p^"^ > 55 GeV/c a difference of about 4% between data and simulation is 
observed. At p^"^ < 55 GeV/c HERWIG differs significantly from the data and PYTHIA in 
the forward region. Due to this large discrepancy we derived the MC corrections from 
the PYTHIA MC. Further discussions of these corrections and HERWIG MC can be found in 

Sec. cni 

6.2 Uncertainties 

If this method was perfect the dijet balance applied to samples after these ?7-dependent 
corrections would be a fiat line at 1. Figures ^1 to 1201 show the data and PYTHIA MC 
after corrections. It is seen that the corrections indeed flatten out the response as desired. 
The remaining discrepancies from a fiat distribution are due to the limitations of the 
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R.g=0.4: • Data a Herwig o Pythia 

Figure 15: Dijet balance, Pdijet = pj^°'"^ / Pt^^^"^^ , as a function of r]jet in data, HERWIG 
and PYTHIA MC samples for Rjet = 0.4 jets. Shown are the corrections for jet-20, jet-50, 
jet-70 and jet-100 jet samples, corresponding to 25 < p^"^ < 55 GeV/c, 55 < p^'^ < 75 
GeV/c, 75 < p^'^ < 105 GeV/c and p^'^ > 105 GeV/c, respectively. The lines show the 
interpolation between the individual measurements used for correcting jets. 



parameterization of the rj- and pt- dependence of the correction and are taken as part of 
the systematic uncertainty of the corrections. 

The systematic uncertainties are determined by varying the event selection require- 
ments and the fitting procedure. Specifically, we varied the cut on the pr of the 3rd jet 
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R.g=0.7: • Data a Herwig o Pythia 

Figure 16: Dijet balance, Pdijet = pj^°'"^ / Pt^^^"^^ , as a function of r]jet in data, HERWIG 
and PYTHIA MC samples for Rjet = 0.7 jets. Shown are the corrections for jet-20, jet-50, 
jet-70 and jet-100 jet samples, corresponding to 25 < p^"^ < 55 GeV/c, 55 < p^'^ < 75 
GeV/c, 75 < p^'^ < 105 GeV/c and p^'^ > 105 GeV/c, respectively. The lines show the 
interpolation between the individual measurements used for correcting jets. 



and the significance of I^t- Any deviation of the dijet balance Pdijet from unity is taken 
as a systematic uncertainty. 

The overall systematic uncertainty from the sources described above are summarized 
in Table El In the low pr region minimum bias data and MC are used up to pt <15 GeV/c 
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R.g=1.0: • Data a Herwig o Pythia 

Figure 17: Dijet balance, Pdijet = pj^°'"^ / p^^^^'^^ , as a function of r] jet in data, HERWIG and 
PYTHIA MC samples for Rj^t = 1-0 jets. Shown are the corrections for the jet-20, jet-50, 
jet- 70 and jet-100 jet samples, corresponding to 25 < pff"^ < 55 GeV/c, 55 < p^'^ < 75 
GeV/c, 75 < p^'^ < 105 GeV/c and p^'^ > 105 GeV/c, respectively. The lines show the 
interpolation between the individual measurements used for correcting jets. 



and single-tower-5 data and the MC jet sample with px > 10 GeV/c for 15 < pr < 25 
GeV/c. Here the systematic uncertainties are largest. 
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Figure 18: Dijet balance, Pdijet = pj^°^^ / p^^^^^^ , as a function of r]jet in data and PYTHIA 
MC samples for cone size Rjet = 0.4 after applying the rj-dependent corrections. Shown are 
the jet-20 (top left), jet-50 (top right), jet- 70 (bottom left) and jet-100 (bottom right) jet 
samples, corresponding to 25 < p^'^ < 55 GeV/c, 55 < p^"^ < 75 GeV/c, 75 < p^""^ < 105 
GeV/c and p'^'^ > 105 GeV/c, respectively. 
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Figure 19: Dijet balance, Pdijet = pj^°^^ / p^^^^^^ , as a function of r]jet in data and PYTHIA 
MC samples for cone size Rjet = 0.7 after applying the rj-dependent corrections. Shown are 
the jet-20 (top left), jet-50 (top right), jet- 70 (bottom left) and jet-100 (bottom right) jet 
samples, corresponding to 25 < p^'^ < 55 GeV/c, 55 < p^"^ < 75 GeV/c, 75 < p^""^ < 105 
GeV/c and p'^'^ > 105 GeV/c, respectively. 
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Figure 20: Dijet balance, Pdijet = pj^°'^^ / p^^^^^^ , as a function of r]jet in data and PYTHIA 
MC samples for cone size Rjet = 1.0 after applying the rj-dependent corrections. Shown are 
the jet-20 (top left), jet-50 (top right), jet- 70 (bottom left) and jet-100 (bottom right) jet 
samples, corresponding to 25 < p^'^ < 55 GeV/c, 55 < p^"^ < 75 GeV/c, 75 < p^""^ < 105 
GeV/c and p'^'^ > 105 GeV/c, respectively. 
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Table 2: Systematic uncertainties of the ?7-dependent corrections versus pj: and \rjjet 



\rj\ range 


0.0 - 


0.2 


0.2 - 


0.6 


0.6 - 


0.9 


0.9- 


1.4 


1.4- 


2.0 


2.0 - 


2.6 


2.6 - 


- 3.6 


PT < 12 GeV/c 


1.5 


% 


0.5 


% 


1.5 


% 


2.5 


% 


1.5 


% 


5.0 


% 


7.5 


% 


12 < pT < 25 GeV/c 


1.5 


% 


0.5 


% 


1.5 


% 


1.5 


% 


1.5 


% 


3.0 


% 


6 


% 


25 < PT < 55 GeV/c 


1.0 


% 


0.5 


% 


1.0 


% 


1.0 


% 


0.5 


% 


1.5 


% 


6 


% 


> 55 GeV/c 


0.5 


% 


0.5 


% 


0.5 


% 


0.5 


% 


0.5 


% 


1.5 


% 


6 


% 
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7 Absolute Jet Energy Scale 



The absolute correction aims to transform the jet energy measured in the calorimeter 
into the energy corresponding to the underlying particle jet. After this correction the 
energy scale of a jet is independent of the CDF detector. Since the calorimeter simulation 
has been optimized to reproduce the measured single particle response, we rely on the 
simulation to derive corrections over a large range of jet transverse momenta. 

The accuracy of this method depends on how well jets are modeled by the simulation. 
In particular it depends on the multiplicity and px spectrum of the particles inside a jet 
and on the response of the calorimeter to an individual particle. These two components are 
tested separately and propagated into a systematic uncertainty on the absolute correction. 

In principle, this correction depends on the initial parton type, e.g. a jet originating 
from a gluon has on average a larger multiplicity than one originating from a quark 
and thus may require a different correction. However, the corrections are not derived 
separately for each parton type since the parton type in the data is a priori unknown. 



7.1 Correction Procedure 

The absolute jet energy is defined as the most probable value for a jet transverse momen- 
tum, pi^^, given a particle jet with a fixed value of pjf''**'^'*^. The corresponding probability 
density function, dP, is parametrized as a function of Ap^ = pPf^*"^'^ _ ^Jf* according to 

/(Apt) = -^=: r[e ' + N2e~^^^^' ] (26) 



^{ai + N2a2) 

One Gaussian function describes the tails while the other one reproduces the core of 
the distribution. Their relative contributions are determined by the normalization of the 
second Gaussian, N2. 

The parameters /xi, ai, fi2, o"2 and N2 depend on p^''^^'^'-'^ as shown by the different 
shapes of the histograms in Fig. |^ The dependence is modeled by a linear parameteri- 
zation of the form: 

^1,2 = <2 + <2Pr*^'^'^ (27) 

/^1,2 = + /X?,2Pr*-'^ (28) 

N2 = + N^p^T^'''"'' (29) 

using a total of 10 parameters. 

The number of jets n with a pt between pj^^ and p^^* + dpj^^ and given particle jet with 

1 ; particle i particle . i particle • • i ji i j.- c jn/ jet particle\ 

a Pt between andpy +dptrp is given by the convolution of dF(P2-j ,PT,i ) 

with the Pt spectrum of the particle jets dV (jf^l^^'^^'^) . That is. 
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/ jet particles I jet i particle / particles in/ i^t particles /onN 

n{prj, ,pj. )dp^ dprp = n[p^ ) x aPyprp ) [SO) 

The likelihood has the form: 

£ = UidVip'Pf''') X dV{plr'\ VtT'^'') (31) 

where the product goes over all particle jets. An unbinned likelihood fit is used to extract 
the parameters of /(Apr) parameters, maximizing the logarithm of the likelihood 

N N 

log C = Y: log dP(/^[*-'^) + log /(Apt.) (32) 

i=l i=l 

where N is the number of particle jets. The first term is independent of the parameters 
to determine and is thus ignored. 

The dijet PYTHIA MC samples, described in Sec. EJ are used to calculate the likelihood. 
Jets are reconstructed at the calorimeter and particle level using the standard CDF jet 
clustering algorithm with cone radii of 0.4, 0.7 and 1.0. Jets are required to be in the 
central region (0.2 < |?7| < 0.6) and to be one of the two leading jets. Each particle jet 



is required to match its closest calorimeter jet within AR = y(A0)^ + (Ar/)^ < 0.1. In 
total we selected about 50,000 particle jets matched to calorimeter jets, with pt ranging 
from to 600 GeV/c. The difference between the particle jet px and the calorimeter jet 
Pt is shown in Fig. |^for four example px ranges. The distributions are not centered at 
zero as expected and have widths that change with pt according to Eq. 1^ 

The absolute correction is shown in Fig. |221 for the three cone radii. At pj^^ = 8 
GeV/c the correction factor is about 1.4 and decreases toward high p^*^* to an asymptotic 
value of about 1.12. At high p^?* the corrections are independent of the cone size while at 
low pg?* a slight dependence is observed. For p^^* < 8 GeV/c, a large fraction of jets are 
not reconstructed since the observed single calorimeter tower energy often falls below the 
1 GeV/c seed tower threshold. In these cases, it is not possible to establish a mapping 
between calorimeter and particle jets and no correction is derived. 



7.2 Uncertainties 

The validity of the probability function used to determine the absolute correction depends 
on how well the MC simulation models the jet response in data. The treatment of the 
jet response as a convolution of the single particle response with the pt spectrum of the 
particles in a jet allows for propagation of the uncertainties of the individual components 
to an uncertainty of the absolute correction. 

Given the calorimeter response R{p) to charged and neutral particles with momen- 
tum p, the average expected response Rave, or jet response, can be calculated from the 
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Figure 21: ApT = p^"^**'^'*^ — p^^^ for Rjf,t=0.4 jets matched using AR < 0.1 for different 



Jet 



particle i • 

Pj, bins. 



generated particles using 



P _ Ef=iPiR{Pi) 

l^i=l Pi 



where is the number of particles inside the jet cone with momenta pi. The uncertainties 
on the absolute corrections is defined as the differences between data and simulation in the 
calorimeter response to single particles, R{p), the differences in the momentum spectrum 
of the particles, p, and from the stability of the calorimeter calibrations in data. They 
are addressed in the following sections. 
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Figure 22: Absolute corrections for different cone sizes as a function of calorimeter jet pt- 
The solid line sliows the corrections for cone size 0.4, the dashed for 0.7 and the dotted 
for 1.0. 



7.2.1 Single Particle Response Simulation 

The measured calorimeter response R{p) = {E/p) (see Sec. E} for hadrons is parameter- 
ized as follows: 

j9 < 20 GeV/c : R{p) = 0.70 + 0.09 tanh(2.13(log(p) - 0.93)) 

p>20 GeV/c : R{p) = 0.70 + 0.14tanh(0.49(log(p) - 1.15)) (34) 

using the data in Fig. [7| For electromagnetic particles (electrons and photons) we set 
R{p) = 1.0. 

The relative uncertainty on the jet energy scale response, AE/E, is then 
AE/E = {E-E^)/E 
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N N 

1=1 1=1 

N 

= TY.PMPi)^(E/p)/E (35) 

i=l 

where A{E/p) is the uncertainty on (E/p) given in Tabled For hadronic particles it is 
between 2.5% and 4%. However, since only 70% of the jet energy arises from hadronic 
particles the uncertainty on the jet energy is only 1.8 — 2.8%. The response function is 
then convoluted with the particle momentum spectra in jets to estimate the systematic 
uncertainty as function of the jet energy. The uncertainty on the simulation of the re- 
sponse of electromagnetic particles is 1.7%. Since about 30% of the jet energy is due to 
electromagnetic particles, this results in a systematic uncertainty of 0.5% on jet energy 
scale, independent of jet pt- The uncertainties are shown versus jet pt in Fig. 123 

7.2.2 Fragmentation 

Uncertainties related to the particle momentum spectrum in a jet originate from the 
modeling of hadronization effects using PYTHIA and HERWIG as well as from the estimate 
of track reconstruction efficiencies in data and detector simulation. The transverse mo- 
mentum spectrum of tracks in data is corrected for inefficiencies as follows. The track 
reconstruction efficiency is measured by embedding simulated tracks inside jets in data 
events, after tuning the simulation of COT hits to distributions observed in data, and 
then parameterized versus jet pt, track momentum and distance of track from the jet 
core. This parameterization is used to correct the data for any inefficiencies |31j . 

The average number of tracks is measured as a function of the track momentum for 
different values of p^^*. Only tracks that are within the jet cone are associated with the 
jet. To account for underlying event contributions, tracks from the region transverse to 
the leading jet are subtracted event by event. Only events with exactly one reconstructed 
2;-vertex are used to reject events with additional pp interactions. Figure |2S1 shows a 
comparison of the track momentum spectra in data with PYTHIA and HERWIG simulation 
for four values of p^*^*. The data are generally in good agreement with the MC, apart from 
some discrepancies at very low track momenta. 

The systematic uncertainty on the jet pt due to the particle multiplicity is calculated 
from the differences between the measured and simulated average calorimeter response. 
Rave for a fixed single particle response R{p). Figure shows the jet response for data, 
PYTHIA and HERWIG. Note that this is an indirect measurement of the energy inferred 
from the track momenta. If the response of the calorimeter, R{p), was one, as is the case 
of electromagnetic particles, this quantity would be unity. The deviation from one thus 
quantifies the fraction of the jet energy measured due to the low calorimeter response to 
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Figure 23: Comparison of the particle momentum spectra in dijet events between data 
with PYTHIA and HERWIG. The CDF data have been corrected for track reconstruction 
efficiency. 



charged hadrons. For a fixed R{p) any difference between data and simulation can only 
arise from a difference in the momentum spectrum. 

The top plot of Fig. 1^ shows Rave versus p^^* for data taken with an trigger threshold 
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Figure 24: Top; Jet response Rave for data (closed circles), PYTHIA (open circles) and 
HERWIG (open triangles) for Rj^t = OA jets as function ofp'r^^. Bottom: Difference between 
data and PYTHIA (open circles) and data and HERWIG (open triangles). 



of 20 GeV/c (see Sec. i)) compared to PYTHIA and HERWIG simulation. At pi^^ = 15 GeV/c 
about 20% of the jet energy is not measured. The response improves with increasing p^*^* as 
expected: e.g., for p^*^* = 120 GeV/c only 14% is not measured. The data response is about 
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0.5% higher than that of the simulation. The bottom plot shows the difference between 
data and simulation. The largest observed difference of 1% is taken as a px independent 
systematic uncertainty on the jet energy scale due to differences in the particle momentum 
spectrum (see Fig. 

7.2.3 Stability of the Calorimeter Energy Scale 

The simulation is tuned using Run 11 data collected during a fixed period of time. The 
calorimeter calibration is kept constant to within 0.5% as described in Sec. 12.31 This 
value is taken to be an additional systematic uncertainty. 

7.3 Summary 

A summary of all the uncertainties is shown in Fig. 1221 It rises from 2% at low p^^* to 
3% at high p^^*. The dominant uncertainty arises from the uncertainty on the simulation 
of the calorimeter response to charged hadrons. The individual uncertainties are added 
in quadrature to give the total uncertainty. The uncertainties apply to all cone sizes. 
Further studies of the validity of the systematic uncertainties are presented in Sec. ^1 
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Figure 25: Systematic uncertainty on the jet energy scale due to the calorimeter calibra- 
tion and simulation. The solid line shows the total uncertainty and the other lines show 
individual contributions. 
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8 Multiple pp Interactions 



At high instantaneous luminosities more than one pp interaction occurs in the same bunch 
crossing at the Tevatron due to the large pp cross section. The number of pp interactions 
per bunch crossing, iV, follows a Poisson distribution with a mean (A^) which depends 
linearly on the instantaneous luminosity. For the Run II configuration of the Tevatron, 
the average number of interactions is about one at L = 0.4 x 10^^ cm~^s^^ and increases 
to three at L = 1 x 10^^ cm"^s"^ and eight at L = 3 x 10^^ cm~^s"^. For the data taken 
up to September 2004 the instantaneous luminosity ranges between 0.1 x 10^^ cm~^s~^ 
and 1 X 10^^ cm"^s"^. 

These extra pp interactions increase the energy of the jets from the hard scatter if 
their final state hadrons accidentally overlap with the jets. This extra energy therefore 
needs to be subtracted from the jet energy. 

8.1 Correction Procedure 

The number of reconstructed 2- vertices, N^t^i is the best estimate of the number of in- 
teractions in a bunch crossing. Vertices are reconstructed using the intersections of the 
tracks with the beam line. Figure 1211 shows the number of vertices as a function of the 
instantaneous luminosity in the first 350 pb~^ of CDF data in — > eVe candidate events. 
The mean number of 2;- vertices, (N^tx), is also shown. At low vertex multiplicity the 
expected linear correlation is observed. 

The efficiency of the vertex finding algorithm depends on the track multiplicity. It is 
about 80% for minimum bias events, 98% for W — >■ ei^e events and greater than 99.9% 
for ti events. These efficiencies have been determined in MC samples and verified using 
the fraction of W ^ eu data events whose vertex is within 5 cm of the ^-position of the 
electron track at the beam line. 

The average transverse energy in a cone is measured using the minimum bias data 
sample. The cone is defined using a seed tower randomly selected in the central calorimeter 
region 0.2 < {rjl < 0.6. The transverse energy, Ej^, in this cone is measured as a function 
of the number of vertices for three cone sizes. Figure EH shows Et versus Nytx- The data 
show the expected linear behavior. 

The data are parameterized using a fitted straight line with coefficients given in Table 
EJ The slope parameters give the extra transverse energy per interaction as a function of 
Nytx- Ideally, the intercept should be zero. The intercept values are in all cases close but 
are consistently larger than zero as a consequence of the vertex finding inefficiency. Note 
that the slope and intercept values measured for the three cone sizes are the same when 
taking into account the cone areas. 
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Figure 26: Left: Instantaneous luminosity for W ^ ev events measured up to September 
2004. Riglit: Mean number of reconstructed vertices in W —>■ eu events versus tie 
instantaneous luminosity. Also sliown is a straiglit line fit. 



Table 3: Intercepts and slopes of the multiple interaction correction for the three cone 
sizes. 



Fit parameter 


Cone 0.4 (GeV) 


Cone 0.7 (GeV) 


Cone 1.0 (GeV) 


intercept 


0.006 ±0.001 


0.018 ±0.002 


0.036 ±0.002 


slope 


0.356 ±0.001 


1.056 ±0.001 


2.153 ±0.002 



8.2 Uncertainty 

The validity of this method depends on two aspects of the vertex finding algorithm: 
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Figure 27: {Et) versus the number of z-vertices for Rjet = 0.4 (top left), Rjet = 0.7 (top 
right) and Rjet = 1-0 (bottom left). A linear fit is also shown. 



• Vertex reconstruction efficiency: The efficiency of finding vertices from ad- 
ditional interactions may depend on the topology of the hard interaction. Any 
inefficiency will result in a steeper slope parameter and a larger intercept. 
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• Vertex fake rate: In events with high occupancy it may happen that fake vertices 
are found, i.e. vertices are reconstructed in ^-positions where no interaction took 
place. This fake rate also depends on the event topology since the probability 
of confusion in both the tracking and the vertex finding increases with increasing 
number of tracks. 

The impact of these effects is tested by repeating the multiple interaction measurement 
using several samples: W — eu^, minimum bias and a jet sample with thresholds of 100 
GeV. Figure OHl shows no indication for any dependence on the instantaneous luminosity 
or on the topology of these samples. However, with the current statistical precision a 15% 
effect cannot be excluded, and it is taken as systematic uncertainty. This value of this 
uncertainty corresponds to 50 MeV for Rjet = 0.4, 150 MeV for Rjet = 0.7 and 300 MeV 
for Rjet = 1.0 per additional interaction. 
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Figure 28: Slope parameter of the multiple interaction correction for Rjf.t=0.7 versus 
instantaneous luminosity in W — »• eu, minimum bias, as well as in a jet sample witli Et 
threshold of 100 GeV. 
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9 Out-of-Cone Energy and Underlying Event 



It is often desirable to reconstruct the energy of the original parton rather than the energy 
of the jet, e.g. for the measurement of the top quark mass or the search for the Higgs 
boson, where parton energies are used to compute the invariant mass of the decaying 
products. 

The reconstruction of the parton energy from the particle jet energy is subject to 
several difficulties. A fraction of the parton energy can be lost from the jet cone due to 
final state gluon radiation (FSR) at large angles with respect to the parent parton or due 
to particles exiting the cone either in the fragmentation process or due to low pt particles 
bending in the magnetic field. This energy is called "Out-of-Cone" (OOC) energy. On the 
other hand the particle jet can also have contributions not related to the actual mother 
parton of the hard interaction of interest defining the jet, such as particles from the initial 
state gluon radiation (ISR), or particles from spectator partons with color connection to 
the other partons of the proton ("beam-beam-remnant", BBR). These two contributions 
are called "Underlying Event" (UE). 

Final state radiation and hadronization effects are correlated with the primary jet 
direction and the jet energy and are expected to decrease with increasing distance from 
the jet core. The UE is thought to be uncorrelated with the direction of the outgoing 
parton and thus independent of the distance from the jet in r] — (j) space and almost 
independent of the jet energy. 

In this section, we derive corrections for the OOC energy and the UE simultaneously 
using PYTHIA dijet MC samples. As in the case of the absolute corrections, the corrections 
are obtained using jets with 0.2 < \r]jet\ < 0.6 since any ?7-dependence of the OOC energy 
is taken into account by the relative corrections. The corrections are solely determined 
from MC simulation at particle generator level independent of the CDF detector. The 
systematic uncertainties of the OOC and UE corrections are derived from comparisons of 
the energy measured in calorimeter towers in certain annuli around the jet cone with the 
simulation based on PYTHIA and HERWIG. 



9.1 Correction Procedure 

The OOC and UE corrections are obtained from PYTHIA dijet samples using particle jets 
which match a primary parton within AR < 0.4. We parameterize the difference of the 
energy between the particle jet and the parton using the same method as for the absolute 
corrections (see Sec. 17. ip . Figure QUI shows _ ^rticie different parton momenta 

and Rjet = 0.4, 0.7 and 1.0. The energy outside the jet cone depends strongly on the cone 
size. The OOC corrections g^j-e shown in Fig. EDI For the smallest cone 

size, Rjet = 0.4, it is about +18% at p^f'^**'^''^ = 20 GeV/c. For the largest cone size, the 
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Figure 29: Normalized distributions of p^^^°" _ ^jmrucie 
0.4 ftoD). 0.7 fmiddle) and 1.0 fbottom). 



for different p^^^"" for cone sizes 



correction is negative: —6% at p^f**'^''^ = 20 GeV/c, corresponding to 1.2 GeV/c. This 
shows that at small cone sizes the OOC losses dominate over the energy increase due to 
the UE, and at large cone sizes the extra energy from the UE is larger than the OOC 
losses. We have estimated that the UE transverse energy is about 0.4 GeV, 1.1 GeV and 
2.2 GeV for cone sizes of 0.4, 0.7 and 1.0, respectively (see Sec. I9.2.2|) . 
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(solid line), 0.7 (dashed line) and 1.0 (dotted line). 



9.2 Uncertainties 
9.2.1 Out-of Cone Energy 

We determine the uncertainty on the OOC energy using 7+jets samples. The reference 
energy scale is the photon pt which serves as an estimator of the corrected jet pt, i.e. 
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The transverse energy around a jet of cone size Rjet is measured by adding the trans- 
verse energy in towers within the annulus defined by radii ri and r2 around the jet axis, 
that is 



PT[ri - rsj 



\ 



N 



N 



(36) 



i=l 



1=1 



where N is the number of towers for which ri > y (?7jet — rjiy + {(j)jet ~ > ^"2- Figure 
1^ shows pri^i — in data, PYTHIA and HERWIG for different jet annuh. The shapes of 
the data and simulation distributions agree rather well, and in the following the mean 
value is used to quantify any disagreement between them. 
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Figure 31: Normalized distributions of the momentum in different annuli outside tlie jet 
cone for data, PYTHIA and HERWIG 7+jets events. 



Figure 1221 shows the difference between data and simulation for the mean values of 
Pri^i — ^2) as a function of p^"^^. The largest difference is observed at low p^^^ and is 
about 4%. The systematic uncertainty is defined as the largest difference between data 
and either PYTHIA or HERWIG, and is parameterized as a function of p^^^ but independent 
of the jet cone size. Since this measurement is made at the calorimeter tower level and 
we apply this correction to a jet after absolute correction, the uncertainties shown in Fig. 
021 are multiplied by the factors 1.1, 1.35 and 1.6, for i?jet=0.4, 0.7 and 1.0, respectively. 



58 



These factors were determined from PYTHIA by comparing the particle and the calorimeter 
energy inside the annuli around the jet cone. 



0.06 
0.04 
0.02 


-0.02 
-0.04 
-0.06 

0.06 
0.04 
0.02 


-0.02 
-0.04 
-0.06 

0.06 
0.04 
0.02 


-0.02 
-0.04 
-0.06 



2t 



2d_ 



Cone 1.3 - Cone 0.4 




_ao 



Cone 1.3 - Cone 0.7 



4n 



_12D CBoiLlGfiMaO 



Cone 1.3 - Cone 1.0 

^^^ ^^^^ 6 



20 



40 



60 



80 



100 



120 
pT' (GeV/c) 



OOC uncertainty a Herwig 



Pythia 



Figure 32: Systematic uncertainty of OOC corrections for different cone sizes. The sys- 
tematic uncertainty is taken as the largest difference between data and either PYTHIA or 
HERWIG and multiphed by the particle-jet /calorimeter- jet factors. 
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By considering alternative generators various modeling uncertainties contributing to 
the systematic error are taken into account. HERWIG and PYTHIA have very different beam- 
beam remnant contributions. Furthermore, they differ in the modeling of QCD radiation 
and fragmentation. 

9.2.2 Underlying Event 

Another source of systematic uncertainty comes from the varying UE in different physics 
processes. To first order these dependencies are taken into account by the MC generators, 
e.g. PYTHIA has been tuned to describe the UE in the data (PYTHIA Tune A, The 
UE uncertainties are derived from comparisons of the UE in data, PYTHIA Tune A and 
HERWIG [n2|- This comparison is done using tracks with px > 0.5 GeV/ c that are separated 
from the leading jet in azimuth by 60° < A$(jet, track) < 120°, which is referred to as the 
"transverse region" . It is mostly sensitive to ISR and multiple parton interactions. Figure 
l33l shows the average momenta of the tracks in the transverse region versus the leading 
jet pt- The data agree well with PYTHIA but differ by up to 30% from HERWIG. This value 
is taken as the relative systematic uncertainty. As a further cross check. Fig. EHl shows 
also the corresponding transverse momentum spectrum simulated by ISAJET which 
has an alternate hadronization model. To get an estimate of the absolute UE uncertainty 
we use the energy measured in minimum bias data for N^tx = 1 as shown in Fig. EZI The 
numbers are 0.4 GeV, 1.1 GeV, and 2.2 GeV for Rj^t = 0.4, 0.7 and 1.0, respectively, 
which translate to UE uncertainties of 0.11 GeV, 0.32 GeV and 0.66 GeV. We have also 
compared the average transverse momenta between data, PYTHIA and HERWIG in 7+jet 
and Z-iet events and find a similar agreement. 

The resulting contribution to the systematic uncertainty of the jet energy scale is 
about 10% at px = 10 GeV/c and decreases to about 2% at px = 70 GeV/c. It can 
be further improved by a more detailed comparison of data and simulation, thus leading 
to a better understanding of the physics effects, or using improved versions of the MC 
generators as JIMMY or PYTHIA 6.3. 

9.2.3 Splash-Out 

The OOC energy refers only to the energy lost outside the jet cone up to i?jet=l-3. In 
PYTHIA MC samples we measured that an additional energy of 0.5 GeV falls outside a 
cone of 1.3. We take half of this energy as systematic uncertainty, i.e. 0.25 GeV, and 
refer to it as "splash-out" uncertainty. 
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Figure 33: The average transverse momentum of charged particles in the transverse region 
as described in the text versus leading jet pt \<^^ calculated using tracks in Rjet=0.7. The 
data are shown as points and are compared to the predictions from PYTHIA (dashed-dotted 
line), HERWIG (solid line) and ISAJET (dashed line). 
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10 Validation of the Jet Energy Scale Determination 



Several consistency checks and further studies are presented in this section. The jet 
energy corrections, which are mostly derived from dijet samples are applied to 7-jet, and 
Z-jet, and tt events to verify the validity of the corrections and systematic uncertainties. 
Furthermore, we present additional studies on the ?7-dependent corrections. 

10.1 Test of the Jet Corrections 
10.1.1 Using 7-jet Events 

The 7-jet data sample is ideal for studying the jet energy scale. The photon energy pj, is 
measured accurately in the CEM calorimeter and thus provides a perfect reference for the 

jet 

jet energy. At tree-level the jet energy should always balance the photon energy: = 1. 
Even in the presence of higher order QCD corrections, which spoil this exact balancing, 
the comparison of the jet energy measurement in these events provides an excellent testing 
ground for the uncertainty on the jet energy scale. In particular, the extent to which the 
data agree with the MC simulation tests the systematic uncertainties in the jet energy 
measurement. 

Photons are selected with > 27 GeV/c and \ri'^\ < 0.9. Jets fragmenting into a tt" 
or rj can decay into photons which constitute a significant background. We require less 
than 1 GeV of extra transverse energy in a cone of radius 0.4 around the photon in the 
calorimeter and less than 2 GeV/c for the scalar sum of the track px values inside the 
cone. We also apply cuts on the shower shape and the number of clusters in the CES 
detector. However, even after these cuts there is still a residual background of about 30% 
at p^ = 27 GeV/c. The background estimate is based on the number of hits in the CPR 
detector jSH]- We estimate the 7-jet balance separately for the signal and background 
and find the results to be consistent for to within 1%. However, the agreement strongly 
depends on the cuts used in the analysis, and for looser cuts we observe differences of up 
to 5%. 

Further cuts are applied to reduce the effects from QCD radiation: 

• The photon and the jet are required to be back-to-back in azimuthal angle: ^(f){'y, jet) > 
3 radians. 

• The event has no more than one jet with pt > 3 GeV/c and |?7| < 2.4. 

Furthermore, only events with one reconstructed vertex are used and thus no correction 
for multiple pp interactions is necessary. 

The 7-jet balance is shown versus jet rj after applying the ?7-dependent corrections 
in Fig. El for jet with Rjet = 0.4. The data, PYTHIA and HERWIG show no residual 
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dependence on rjjet as desired. Note that the pt balance between the jet and the photon 
is not expected to be zero at this stage of the correction procedure. The overall scale 
difference between data, PYTHIA and HERWIG will be discussed in Sec. I1U.2I 
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Figure 34: px balance, ^ 1, in data (full circles), PYTHIA (open circles) and HERWIG 

(open triangles) as function of rjjet for Rj^t = 0.4 



Next, we apply also the absolute correction to the jet momenta making the jet pr 
independent of the calorimeter response. Figure shows the resulting pt balance for 
data, PYTHIA and HERWIG for Rj(,t=OA. For comparison, also the pr balance calculated 
using particles at generator level without detector simulation is overlaid. Tables |3] and El 
summarize the mean and the width obtained from fits of Gaussians to these distributions 
within the range -0.4 to 0.4 for all jet cone sizes. For jets from data with Rjet = 0.4 the 
mean is measured to be about 9% lower than the photon. Compared to the data, the 
mean found in PYTHIA is about 2% higher and for HERWIG it is 2% lower. Generally, the 
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Pt of the jet is smaller than the pr of the photon due to the energy lost outside the cone. 
In fact, one observes that for the larger cone sizes, in particular i?jef = 1-0, the mean 
is much closer to since there is nearly no energy lost outside the cone. For particle 
jets, the mean values of PYTHIA and HERWIG agree with the respective calorimeter jets to 
within 1%, which proves the validity of the absolute correction procedure. The observed 
differences between the generators reflect the different modeling of the underlying physics 
process. This difference is largest for i?jei = 0.4. However, the data generally lie between 
PYTHIA and HERWIG and agree to within 2% with both. 



Table 4: Mean value of Pr^ /p'r ~ 1 after //-dependent and absolute energy correction, for 
data, PYTHIA, and HERWIG for i?jet=0.4, 0.7 and 1.0. For PYTHIA and HERWIG, the values 
are given also for particle jets. 



Sample 




Rjet—^-'^ 




Calorimeter jets 


Data 

PYTHIA 

HERWIG 


-0.088 ±0.001 
-0.070 ±0.001 
-0.108 ±0.001 


-0.016 ±0.001 
-0.015 ±0.001 
-0.043 ±0.001 


0.022 ±0.001 
-0.002 ±0.001 
-0.024 ±0.001 


Particle jets 


PYTHIA 
HERWIG 


-0.078 ±0.001 
-0.113 ±0.002 


-0.037 ±0.001 
-0.061 ±0.002 


-0.009 ±0.001 
-0.019 ±0.002 



From Table El we note that the data resolution is around 4-7% worse than HERWIG and 
12% worse than PYTHIA. We observe that HERWIG has a wider resolution than PYTHIA for 
both calorimeter and for particle jets. 

After applying all corrections (?7-dependent , absolute, OOC±UE) we obtain the 7-jet 
balance as shown in Fig. EHl Table IHl contains the corresponding mean values derived 
using a fit of a Gaussian to data and MC distributions for all three cone sizes. Data and 
MC agree with zero to within 2% except for HERWIG for a cone size of i?jet=0.4. The 
differences between data and simulation are equal to those observed in Table IH since the 
OOC and UE correction were derived from PYTHIA and uniformly applied to all samples. 
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Figure 35: 7-jet balance in data, PYTHIA and HERWIG for Rjet=0.4. Overlaid is the cor- 
responding •y-jet balance on particle level jets (triangles) calculated using particles at 
generator level without detector simulation. The distributions are normalized to 1. 



65 



Table 5: Width of /Pr — 1 after the ?7-dependent and absohite energy correction, for 
data, PYTHIA and HERWIG for i?jei=0.4, 0.7 and 1.0. For PYTHIA and HERWIG, the values 
are given also for particle jets. 



Sample 




-Rjei— 0.7 


-Rjet — 1-0 


Calorimeter jets 


Data 

PYTHIA 

HERWIG 


0.199 ±0.001 
0.176 ±0.001 
0.192 ±0.001 


0.191 ±0.001 
0.171 ±0.001 
0.181 ±0.001 


0.191 ±0.001 
0.169 ±0.001 
0.178 ±0.001 


Particle jets 


PYTHIA 

HERWIG 


0.105 ±0.001 
0.127 ±0.002 


0.095 ±0.001 
0.116 ±0.002 


0.090 ±0.001 
0.111 ±0.002 



Table 6: Mean value of Pr^ /pr — 1 after all corrections, including the out-of-cone energy 
correction for data, PYTHIA, and HERWIG for jet cones of Rj^t = 0.4, 0.7 and 1.0. 



Sample 


i?jei=0.4 


-Rjet=0.7 




Data 

PYTHIA 

HERWIG 


-0.019 ±0.001 
-0.001 ±0.001 
-0.040 ±0.001 


0.010 ±0.001 
0.011 ±0.001 
-0.018 ±0.001 


0.024 ±0.001 
0.000 ±0.001 
-0.023 ±0.001 



10.1.2 Using Z-jet Events 

Another excellent calibration sample are Z events where the px of the Z boson 

provides a reference scale for the jet. The advantage compared to the 7-jet sample is that 
it is nearly free from background contamination, at the expense of smaller statistics. For 
this study we require the jet and the Z boson to be back-to-back, A(^(jet, Z) > 3 radians, 
and no extra jets with > 3 GeV/c and |?7| < 2.4. 

In Fig. inZlwe compare the Z-]ei balance, Pr^ /p^ — 1 in Z ^ e^e~ and Z 
events for data, PYTHIA and HERWIG after all corrections. 

The mean values derived from fits of Gaussians to the distributions between —0.3 and 
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Figure 36: 7-jet balance in data, PYTHIA and HERWIG using Rj(.t=0.4, 0.7 and 1.0 after 
r]-dependent, absolute and OOC+UE corrections. 



+0.3 are given in Table They are reasonably close to for all cone sizes, and the data 
agree with the MC to within the statistical uncertainties of 1%. 
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Figure 37: Z-jet balance in data (closed black circles), PYTHIA (solid red line) and HERWIG 
(dashed blue line) for Rj^t = 0.4 (top left), Rj^t = 0.7 (top right) and Rj^t = 1-0 (bottom 
left) after all corrections. 
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Table 7: Measured mean value of Vpf^ ^ after all corrections, including the out-of-cone 
energy correction for data and PYTHIA calorimeter jets. The result is given for Rjet = 0.4, 
0.7 and 1.0. 



Sample 


Rjet =0.4 


Rjet —0.7 


Rjet —1.0 


Data 

PYTHIA 

HERWIG 


-0.026 ±0.009 
-0.016 ±0.003 
-0.032 ±0.003 


0.007 ±0.009 
0.019 ±0.003 
-0.011 ±0.002 


0.013 ±0.009 
0.015 ±0.003 
-0.009 ±0.003 



10.1.3 Using Dijet Events 

The ?7-dependent corrections have been obtained from the PYTHIA dijet and jet data 
samples. As explained in Sec. IHl the data and MC have different ?7-dependent corrections. 
This section describes the choice of the dijet balancing technique as the r^-dependent 
correction method and the use of only PYTHIA MC to correct all MC samples. 

A different approach to address the ?7-dependence of jet response is the so-called "Miss- 
ing Et Projection Fraction" (MPF) method. This approach was used in Run I by the 
CDF jSHj and DZero [lU] experiments. The MPF is defined as 



rf probe 

MPF = . (37) 

where the vector of the missing transverse energy, I^t, is used to quantify the difference 

1 , probe 1 triqqer ,1 ,1 • probe triqqer • i • j-i j-- j. i_ 1 

between andp^^ rather than usmg —Pr as is done m the dijet balancmg 
method (see Sec. Ej). In an ideal dijet production process with no gluon radiation and 
fragmentation effects, the MPF and the dijet balance methods are equivalent. That is, 

^^^^^^ = 2 ± {MPF) - 2^ = 

However, due to QCD radiation and out-of-cone energy losses and underlying event con- 
tributions this does not exactly hold. In contrast to the dijet balancing method, the 
MPF method does not correct for OOC energy. The reason is that I^t used in the MPF 
method is only affected by energy mismeasurement and has no sensitivity to the energy 
flow between inside and outside the jet cone. On the other hand, the dijet balancing 
method is sensitive to the energy inside the jet cone and will thus implicitly correct for 
an ?7-dependence of the OOC and UE effects. Since we estimate these corrections and 
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systematic uncertainties only in the central region, we choose to use the dijet balancing 
method as the primary correction method. 

However, we use the MPF method to further investigate the discrepancy between 
PYTHIA and HERWIG in Figures IT^HTTI in Sec. El The HERWIG measurements are systemat- 
ically higher than PYTHIA and data by about 10% at 25 < p^" < 55 GeV/c and |?7| > 0.6 
and agree very well at higher GeV/c. The discrepancy is larger for i?jet=0.4 jets than 
for Rjf.t=0.7 and 1.0. To shed more light on the origin of the discrepancy we compare the 
ratio (3{MPF) / (3dijet as a function of r/jet and for two different values oi p'^'^ in Fig. |2H1 
In the central region (3{MPF) / ^dijet is consistent with unity at both values oip^'^ . In the 
forward region [3{MPF) increases with respect to I3dijet- For 25 < pff*^ < 55 GeV/c the 
data are well modeled by PYTHIA MC but large discrepancies are observed in comparison 
to HERWIG MC. At high p'^'' > 105 GeV/c the data are in good agreement with both MC 
generators. Since the ratio P{MPF) / Pdijet is largely independent of the CDF jet energy 
scale we conclude that the observed disagreement is not due to any residual problems 
of the CDF simulation but must originate from a difference in the underlying physics 
between HERWIG, PYTHIA and data for low Pt dijet production. 
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Figure 38: l3{MPF) / Pdijet as a function of rjj^i for Rjet = 0.4 jets in data (circles). 



PYTHIA (upward triangles) and HERWIG (downward triangles) for 25 < p^' 
and p^^'^ > 105 GeV/c. 



< 55 GeV/c 



We also compare the (3dijet for particle-jets with Rjet = 0.4 between PYTHIA and HERWIG 
in Figure inni for pT > 18 GeV/c and find that for PYTHIA it is independent of rjjet while 
for HERWIG it rises with increasing rjjef 

Since this behavior is only found in the dijet samples, we do not consider HERWIG 
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dijet samples for the determination of the ?7-dependent corrections or their systematic 
uncertainties. In 7-jet (see Fig. 2' -jet or tt events no such problems are seen. At this 
moment we do not have any explanation for the differences. It could be due to initial of 
final state radiation, due to the underlying event modeling or many other effects, and it 
will be studied again in future versions of the generators. 
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Figure 39: Pdijet for particle jets as a function of jet rijet for PYTHIA (open circles) and 
HERWIG (open triangles) for pr > 18 GeV/c and for Rjet = 0.4. 



10.1.4 Using W ^ jj Decays in tt Events 

The jet energy scale can be studied using the hadronic decay of resonances with well- 
measured masses such as the W and Z bosons. Unfortunately, the decays of W and Z 
bosons to jets are swamped by multijets QCD background in hadron collider environ- 
ments. One solution is to study hadronic W boson decays inside tt events which have 
relatively small background contamination. This section summarizes the application of 
this technique. For complete details, see [HT] . 

At the Tevatron, top quarks are produced primarily as top pairs and decay to W 
bosons and b quarks nearly 100% of the time. The W bosons in turn decay into lepton- 
neutrino (lu) or quark pairs {qq'). This measurement uses the "lepton+jets" channel of 
tt candidates in which only one of two W bosons decays to lu while the other decays to 
quark pairs. The lepton-|-jets events are selected by requiring one well-identified electron 
or muon, large (-^t) due to the neutrino from the W decay and at least four jets in 
the final state. The missing transverse energy, I^t is measured by the imbalance in 
the calorimeter transverse energy and is required to be greater than 20 GeV. Jets are 
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reconstructed with a radius Rjet = 0.4. The sample is divided into four subsamples with 
various sensitivities for better performance. First, the events are separated based on the 
number of jets that are 6-tagged in the event. The SECVTX [3H] algorithm based on the 
identification of secondary vertices inside jets is used to tag fe-jets. Events with 2-,l- and 
0-tag are considered separately. Furthermore, events with 1-tag are separated based on 
the fourth jet Et threshold. Events in the l-tag(T) category have 4 jets with Et > 15 
GeV, while events in the l-tag(L) category have 3 jets with E^ > 15 GeV and the 4th jet 
with 8 < < 15 GeV. 

Before reconstructing the invariant mass of hadronically decaying W bosons {tTLjj), we 
apply the ?7-dependent, and absolute corrections to jet energies. In addition, corrections 
specific to light quark jets from W boson decays in ti events are applied. To reconstruct 
rrijj, one has to know which of the jets in the final state comes from the W boson decay. 
This problem is dealt with by considering all the dijet combinations that can be made using 
the jets that are not fe-tagged. Only the four highest Et jets are considered. Consequently, 
there can be more than one mass per event that are considered. There are in fact 1,3, and 
6 nijj per event for the 2-tag, 1-tag and 0-tag subsamples, respectively. The distribution 
of rrijj for HERWIG ti events is shown in Fig. El for each event category (with a top quark 
mass (Mtop) of 178 GeV/c^). The mass resolution improves with the number of 6-tagged 
jets present in the event. 

Distributions of constructed from HERWIG tt Monte Carlo with Mtop = 178 GeV/c^ 

(corresponding to the central value of the Tevatron Run I average) with jet energy scale 
values ranging from -3 to +3(Jc, where cXc is the total jet energy scale uncertainty defined 
in Sec. ^2 of this document. Smooth probability density functions are obtained by fit- 
ting the mass distributions as a function of Mtop and jet energy scale using an analytical 
function. Figure 1^ shows the rrijj distribution for various jet energy scale values for the 
2-tag subsample with the fitted templates overlaid. Templates for background events are 
obtained from ly+jets, QCD multijets and single-top MC events. 

The fitted jet energy scale is obtained by comparing the reconstructed mass distribu- 
tions obtained in the data with the signal and background templates using an unbinned 
likelihood fit. The data used in this measurement corresponds to the W — eu^ and 
W — >■ fif^ but additionally requiring at least 4 jets and I^t greater than 20 GeV. The 
systematic uncertainties in this measurement arises from the MC modeling of signal and 
background events that we use to create the templates and thus extracting the jet energy 
scale. We consider uncertainties in the top quark mass (± 5 GeV/c^), amount of initial 
and final state gluon radiation, parton distribution functions, background mass shape and 
general MC modeling. The total systematic uncertainties corresponds to 0.68 Uc. 

This measurement is performed using 318 pb~^ of data that results in a total of 
165 events in the lepton-|-jets sample. The application of the likelihood fit to the data 
yields —0.76 ± 1.00 (stat.) cTc- This means that the data and simulation of ti events is 
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Figure 40: Reconstructed hadronic W boson mass from HERWIG tt events with M^p = 
178 GeV/c^ for 2-tag events (upper-left), l-tag(T) events (upper-right), l-tag(L) events 
(bottom-left) and 0-tag events (bottom-right). The yellow (outer) histograms show the 
mass distributions for all combinations and the blue (inner) histograms show the distri- 
butions only for the correct jet-parton assignments. 



in agreement in one ac- By adding the systematic uncertainties, this result changes to 
—0.76 ± 1.27 cTc- The rrijj distributions reconstructed in the data are shown in Fig. 
The shape of the signal and background MC templates corresponding to the best fit are 
overlaid on top of the histograms. We conclude that the average jet energy scale as 
determined hj W ^ jj decays is in good agreement with the nominal jet energy scale of 
the CDF MC simulation. 

The calibration of the jet energy scale with W — >■ jj decays has been used to measure 
precisely the top quark mass in Run II j^H]- We note that jet energy scale uncertainties 
obtained from this technique are mostly statistical and will improve as more data is 
accumulated. 
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Figure 41: Four rrijj signal templates are shown for jet energy scale (JES) values ranging 
from -ScTc to +3o"c. Overlaid are the fitted parameterizations at each value of jet energy 
scale. 

10.2 Test of the Uncertainties 

We test whether the agreement of data and MC is within the calculated uncertainties for 
all pt and rj bins. Figure 1^ shows the difference of the mean values of the 7-jet balance 
between data and PYTHIA as a function of pt and for six regions of pseudo-rapidity. 
Overlaid is the total systematic uncertainty on the jet energy scale. It is seen that the 
data are modeled well by the simulation at all r] and pt, and that any differences are 
covered by the quoted uncertainties. In Fig. |^ the same comparison has been made 
with HERWIG, leading to the same conclusion. The 7-jet balance for data, PYTHIA and 
HERWIG is independent of p^ and rjjet after applying all the corrections. 

10.3 Summary 

We have shown that the corrections and the systematic uncertainties are valid for several 
control samples. We have found that the transverse momentum of the jet, after all 
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m.. (GeV/c) 

Figure 42: Data nijj distributions for the 2-tag (upper-left), l-tag(T) (upper-right), 
l-tag(L) (lower-left) and 0-tag (lower-right) subsamples. The signal and background tem- 
plate shapes corresponding to the best fit of the jet energy scale cross-check are overlaid 
on the histograms. The value of Mtop has been constrained to 178 GeV/c^. 



corrections, is in balance with the transverse momentum of the 7 and the Z boson in both 
measured and simulated 7-l-jets and Z+jets samples, respectively. We also determined 
the average jet energy scale using W — > jj decays in ti events and found good agreement 
between data and simulation. 

We have also investigated differences between the PYTHIA and HERWIG MC genera- 
tors. With the dijet-balancing technique we observe rather large differences in the plug 
calorimeter region which are also seen for particle jets. For the dijet process, the data do 
not support the behavior of HERWIG. This problem is only observed in the dijet process, 
e.g. not in the 7-jet process. 
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Figure 43: Difference of the 'j-jet balance between data and PYTHIA as a function of pj, 
in six reeions of n^^t. All tliree cone sizes are sliown: 0.4 (blue sauares). 0.7 (red oven 
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Figure 44: Difference of the 'j-jet balance between data and HERWIG as a function of pj, 
in six reeions of r?,»t. All tliree cone sizes are sliown: 0.4 (blue sauares). 0.7 (red oven 



11 Summary of Systematic Uncertainties 



We have presented the systematic uncertainties associated with the jet energy response. 
The systematic uncertainties are largely independent of the correction applied and mostly 
arise from the modeling of jets by the MC simulation and by the knowledge of the response 
to single particles. 

Figure EHl shows the individual systematic uncertainties as a function of jet px in 
the central region, of the calorimeter, 0.2 < \iri\ < 0.6, of the calorimeter. They are 
independent and thus added in quadrature to derive the total uncertainty. 
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Figure 45: Systematic uncertainties as a function of the corrected jet pt in 0.2< \ri\ <0.6. 



For Pt > 60 GeV/c the largest contribution arises from the absolute jet energy scale 
which is limited by the uncertainty of the calorimeter response to charged hadrons. A 
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further reduction of the systematic uncertainties can be achieved by improving the tun- 
ing of the simulation, and by including in situ single track data which recently became 
available, replacing test beam data used so far in the momentum region 7-20 GeV/c and 
probably beyond. 

At low pt the largest uncertainty arises from the out-of-cone energy which can be 
improved by further studying differences between the data and the predictions of PYTHIA 
and HERWIG, and by optimizing the fragmentation and underlying event model of both 
generators. 

Additional uncertainties arise from the fragmentation models, the stability of the 
calorimeter calibration and the underlying event modeling. 
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12 Conclusions 



We have determined a set of corrections to estimate the parton energy from the jet 
energy measured in the Colhder Detector at Fermilab. The cahbration is based on data 
taken between 2001 and 2004 at the Tevatron pp colhder, corresponding to an integrated 
luminosity of about 350 pb^^, and on test beam data. 

These corrections involve several steps and for each step a systematic uncertainty is 
determined. Both the central and forward components of the calorimeter are calibrated 
using test beam and in situ data. The response of jets in the forward calorimeter is 
calibrated with respect to that of the central calorimeter. The shower simulation is in 
particular tuned in detail to the data in the central rapidity region. Using several MC 
generators it has been shown to provide a good description of the energy response of 
various physics processes. The MC simulation is used to derive a correction for the 
calorimeter jet energy response in the central region. Further corrections are made for 
multiple pp interactions, the underlying event and the fractional energy of the parton 
that is not contained within the jet cone. Finally, we have verified that the corrected jet 
energy is a good measure of the initial parton energy using prompt photon and Z events 
and have shown that the various MC generators provide a good description of the data 
within the quoted systematic uncertainties. 

The total systematic uncertainty on the jet energy scale varies between 8% at low jet pt 
and 3% at high jet pr- The dominant sources of systematic uncertainty are the uncertainty 
on the test beam measurements at high energy and the uncertainty in modeling the energy 
flow around the jet cone. 
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